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ABSTRACT 

We show that the spectral and radial distribution of the nonthermal emission of massive, M > 
10 Mq, galaxy clusters may be approximately described by simple analytic expressions, which 
depend on the cluster thermal X-ray properties and on two model parameter, f3 cmc and rj e . /3 core is 
the ratio of the cosmic-ray (CR) energy density (within a logarithmic CR energy interval) and the 
thermal energy density at the cluster core, and r] e (p) is the fraction of the thermal energy generated in 
strong collisionless shocks, which is deposited in CR electrons (protons). Using a simple analytic model 
for the evolution of intra-cluster medium CRs, which are produced by accretion shocks, we find that 
/3 cole — r] p /200, nearly independent of cluster mass and with a scatter A In (3 corc ~ 1 between clusters 
of given mass. We show that the hard X-ray (HXR) and 7-ray luminosities produced by inverse 
Compton scattering of CMB photons by electrons accelerated in accretion shocks (primary electrons) 
exceed the luminosities produced by secondary particles (generated in hadronic interactions within the 
cluster) by factors ~ 500(77 e /?7p)(r/10keV)- 1 / 2 and ~ lSO^e/^XT/lOkeV)" 1 / 2 respectively, where 
T is the cluster temperature. Secondary particle emission may dominate at the radio and very high 
energy (> 1 TeV) 7-ray bands. Our model predicts, in contrast with some earlier work, that the 
HXR and 7-ray emission from clusters of galaxies are extended, since the emission is dominated at 
these energies by primary (rather than by secondary) electrons. Our predictions are consistent with 
the observed nonthermal emission of the Coma cluster for rj p ~ r\ e ~ 0.1. The implications of our 
predictions to future HXR observations (e.g. by NuStar, Simbol-X) and to (space/ground based) 7- 
ray observations (e.g. by Fermi, HESS, MAGIC, VERITAS) are discussed. In particular, we identify 
the clusters which are the best candidates for detection in 7-rays. Finally, we show that our model's 
results agree with results of detailed numerical calculations, and that discrepancies between the results 
of various numerical simulations (and between such results and our model) are due to inaccuracies in 
the numerical calculations. 

Subject headings: acceleration of particles - galaxies: clusters: general - radiation mechanisms: non- 
thermal - X-rays: general 



1. INTRODUCTION 

A nonthermal emission component is observed in sev- 
eral clusters of galaxie s. In most ca ses nonthermal radio 
emission is observed (jFerettil [20051 ). and in some cases 
a nonthermal hard X-ray (H XR) emission compon ent is 
also observed (for review, see iRephaeli et aLll2008l l. The 
radio emission is interpreted as synchrotron radiation, 
thereby suggesting that relativistic electrons and mag- 
netic fields are present in the intracluster medium (ICM) . 

Several models for the synchrotron emission have been 
presented in the literature. These models differ in the as- 
sumptions regarding the origin of the emitting electrons. 
In some models the emitting electrons are secondary elec- 
trons and positrons that were generated by p-p interac- 
tions of a cosmic ra y (CR) proton population w ith the 
ICM (e.g. iDennisor] 119801 : iBlasi fc Colafrancescol 11999( 1 . 
while in others the high energy of the emitting elec- 
trons is acheived by turbulent acceleration of a preexist- 
ing ICM population of no nthermal " seed elect r ons" (sec- 
onda r y or otherwise, e.g. iBrunetti et alj|2001t iPetrosianl 
20011: iBrunetti fc Blasil l2005t ICassano fc Brunettil 120051: 
Cassano et al.1 l2007t IBrunetti et all 12008( 1. Various 



sources of CR protons and electrons were conside red in 
the literature, including active galactic nucl ei (e.g. iKatzl 
119761: iFabian et all 1 19761 iFuTita et all 12007( 1. dark mat- 
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ter bow shocks (e.g. iBvkov et al.l I2000D. ram-pressure 
stripping of infalling galaxies (e.g. Ide Plaa et al.l f2006h 
and shock waves associated with the process of large 
scale s tructure (LSS) formation Ce.g. ILoeb fc Waxmanl 
20001: iFuiita et all l200l iBerrington fc Permed 120031: 



Gabici fc Blasil 120031 : IBrunetti et all l2004t llnoue et alJ 



2005). In this paper we consider CRs produced by LSS 
shocks, and derive predictions for the nonthermal clus- 
ter emission they generate. Detailed analysis, based on 
the predictions of this paper, of the nonthermal radio 
and HXR emission from a sample of clusters suggests 
that CR accel eration in LSS shocks is the main source 
of clu ster CRs (jKushnir et alJl2009b iKushnir fc Waxmanl 

Hoioji . 

The following point regarding the general applicability 
of our results should be emphasized here. The simple an- 
alytic expressions we derive (in §[2} give the non-thermal 
emission from secondary particles (produced by CR in- 
teractions at the cluster core) as a function of /3 COTe and 
of the observed thermal X-ray properties of the cluster. 
These expressions are valid, for a given value of /3 corc , 
regardless of the ori gin of the CRs residin g at the cluster 
core. As shown in (jKushnir et al.l 1200 9). cluster obser- 
vations already allow one to estimate both the average 
value of /3 C orc and its scatter, Aln/3 corc . The values de- 
rived from observations are naturally explained within 
the frame of the model presented here, where cluster CRs 
are produced by LSS shocks, and difficult to explain in 
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other models. Moreover, since alternative models for the 
origin of cluster CRs must reproduce the observationally 
inferred values of /3 core and Aln/3 core , the non-thermal 
secondary emission predicted by such models would be 
similar to that derived here. 

There are two primary populations of L SS shocks, 
which differ in their Mach numbers fsee e.g. | Ryu et al.l 
[200l iPfrommer et all 12001 iSkillman et all 12008ft : High 
Mach number shocks generated by the accretion of gas 
onto filaments and clusters (accretion shocks), and low 
Mach number shocks produced by halo mergers (merger 
shocks). The hypothesis that accretion shocks gener- 
ate high energy CRs is supported by their resemblance 
to collisionless no n-relativistic shock s in the interstellar 
medium (see e.g. iKeshet et al.l [2003). which are gener- 
ally known to ac celerate a power-law distrib ution of high 
energy particles (jBlandford fc Eichler|[l987l) . 

Accretion shocks are expected to accelerate electrons 
to TeV energies. Such electrons, to which we refer as 
primary electrons, lose their energy by inverse Compton 
(IC) scattering of CMB photons on a time scale which 
is much shorter then the cluster dynamical time, leading 
to HXR a nd 7-ray emissi o n at t h e vicinity of the shocks 
see e.g. iDar fc Shavivl 119951: IColafrancesco fc Blasil 



19981: iLoeb fe Waxmanl 120001: IWaxman fc Loebl 120001: 
Totani fc Kitavamal 120001: Kawasaki fc Totanil 12002) . 
Accretion shocks are also expected to accelerate protons 
to high energies. The protons, which do not lose their 
energy on a cluster dynamical time, may be coupled to 
the thermal plasma by magnetic fields and be confined 
within the cluster volume over cosmological times. In- 
elastic collisions of these CR protons with thermal ICM 
protons would produce high energy 7-rays through the 
decay of neutral pions, and high energy secondary elec- 
trons iand__pj3sitrons through the decay of charged pi- 
ons ([Dennisonl 1980; V olk et all 119961 : [Berczinsk v et al.l 
|1997|) . The secondary nonthermal emission is dominated 
by the cluster's core, where the ICM and CR density is 
highest. 

To date, nonthermal emission from accretion shocks 
has not yet been reliably observed. Furthermore, no 
clust er of galaxies has so far been firmly detected in 7- 
rays (|Reimer et~aT1[200l_ or in very high energy (VHE, 
> 1 TeV) 7-ravs (jPerkins et al.l [20061: iDomainko et al.1 
l2007t lPerkinsll2008D (see, how ever, the stacking analysis 
of IScharf fc Mukheried [20021 which revealed a low sig- 
nificance excess of 7-ray emission associated with Abell 
clusters). As we discuss in § [5J future HXR observa- 
tions (e.g. NuStar, Simbol-X) and (space/ground based) 
7-ray observations (e.g. Fermi, HESS, MAGIC, VERI- 
TAS) are expected to provide unambiguous detection of 
nonthermal emission from accretion shocks. 

In order to determine the nonthermal relativistic parti- 
cle population in clusters, one has to follow the hydrody- 
namic history of the gas confined within these structures. 
Nonthermal relativistic particles (CRs) are produced in 
accretion and merger shocks, and the CR energy density 
is modified by adiabatic expansion/compression of the 
gas, and by energy loss due to emission of radiation and 
due to inelastic nuclear collisions. The CR population 
in clusters may also be affected by diffusion. Diffusion 
of CRs on scales ~ 100 kpc would reduce their density 
at the cores of clusters and hence the secondary emis- 
sion. The effects of diffusion are difficult to estimate 



theoretically, mainly due to our ignorance regarding the 
structure of magnetic fields within the ICM. However, 
cluster radio observations imply that diffusion does not 
significantly affect the secondary population, and that 
the diffusion time of 100 GeV CRs over scales > 100 kpc 
within the ICM is not short compared to the Hubble time 
(|Kushnir et al J 120091) . We therefore neglect diffusion in 
the present analysis. 

The evolution of CRs in large-scale structures 
was investigated numerically in several earlier studies 
(IKeshet et al.ll2003UMiniatill2003l: IPfrommer et alll2007l 
2008lHfrE mmc r2008). All studies identified cluster cores 
and cluster accretion shocks as the main sources of non- 
thermal radiation. There are discrepancies between the 
results of some of the studies, regarding the particles 
(primaries or secondaries) that dominate the emission at 
different energies and at different spatial locations. Due 
to the complicated nature of numerical simulations it is 
difficult to trace the origin of the discrepancies, which 
may be largely due to the different values adopted in the 
simulations for the acceleration efficiencies in shocks. We 
will resolve here some of these discrepancies by compar- 
ing the numerical results with the results of our analytic 
analysis. 

Several earlier studies attempted to provide an 
analytic derivation of the nonth e rmal luminosity 
of galaxy clusters (jSaraz ml [19991: iFuiitafc barazm 



2001; Bcrrington & Dernier 



2003 



Araudo et all 120081 : iMurase et al.l [2008|). Thesestudies 



Gabici & Blasi 2002 



did not include a model for the evolution of the space 
and energy distributions of primary and secondary CRs 
due to both accretion/merger shocks and adiabatic ex- 
pansion/compression. Moreover, they did not consider 
IC emission by primary electrons, which, as we show be- 
low, dominates the HXR and 7-ray emission. 

We derive in this paper a simple analytic model de- 
scribing the spatial and spectral distribution of the non- 
thermal emission produced by cluster CRs. We focus on 
massive clusters, M > 10 145 M Q , since their nonthermal 
emission is strongest and may be detectable by current 
and upcoming experiments. We assume that accretion 
and merger shocks produce power-law momentum distri- 
butions of CRs, dlnn/dliip = —a, with a determined 
by the shock Mach numbe r, M, a = 4/(1 — M~ 2 ) — 2 
(jBlandford fc Eichlenfl987l ). and that strong (M 2 > 1) 
shocks deposit a fraction ?7p( e ) of the thermal energy they 
generate in CR protons (electrons). Our description of 
the accretion shocks, which dominate the CR generation 
(see § Ej) , is highly simplified. We assume that the clus- 
ters are spherical, that gas is accreted onto clusters in 
spherically symmetric flows at a rate of ~ M/tn, where 
tn is the (instantaneous) Hubble time, and that accre- 
tion shocks are strong (M 2 ^> 1). We show (in §[3| that 
the results of numerical simulations are consistent with 
an average accretion rate (at low z) of ~ 0.5M/i#. 

The nonthermal emission is dominated by primary 
electrons near the accretion shocks and by secondaries 
at the cluster core. It is determined therefore mainly 
by two parameters, r\ e and /3 corc , the ratio of the CR 
energy density (within a logarithmic CR energy inter- 
val) and the thermal energy density at the cluster core. 
We first derive in § [2] simple analytic expressions for the 
nonthermal emission as function of rj e and /3 C oro, and 
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of the thermal X-ray emission properties of the cluster 
(the charged secondary, e^, emission, which may domi- 
nate at radio wavelengths, depends also on the strength 
of the ICM magnetic field). We also show, based on 
simple physical arguments and crude approximations, 
that /3 CO re — 77p/100 should be expected in massive clus - 
ters (see also lPfrommer et al.l2007tlJubelgas et al.ll2008D . 
The results of our simple model are shown to be consis- 
tent with those of detailed numerical simulations in § [31 
with the exception of deviations which are due to inac- 
curacies of the numerical calculations. In § [4] we com- 
pare our model predictions with the observed nonther- 
mal emission of the Coma cluster. The extended HXR 
emission of Coma is shown to be consistent with our 
model predictions, and the HXR and radio luminosities 
are shown to be consistent with our model predictions 
for i] p ~ rj e ~ 0.1. 

Some limitations of our simplified model should be 
highlighted. Our model does not capture temporal fluc- 
tuations in the mass accretion rate and deviations of the 
accretion flow from spherical symmetry. Since the cool- 
ing time of primary electrons is short compared to the 
cluster dynamical time, these effects will lead to tempo- 
ral fluctuations and to deviations from spherical symme- 
try of the HXR and gamma-ray emission. Our model 
approximately describes therefore only the temporally 
and azimuthally averaged emission of primary electrons. 
Dedicated numerical simulations that constrain the tem- 
poral and spatial fluctuations of the accretion rate will be 
useful for analyzing upcoming observations (as discussed 
in detail in § [5]). The emission produced by secondaries, 
on the other hand, is not sensitive to the spatial and tem- 
poral fluctuations in the accretion of gas. This is due to 
the fact that the cooling time of CR protons, which pro- 
duce the secondaries, is long (of the order of the Hubble 
time). They therefore accumulate in the cluster over its 
age, and their density at the cluster core, which domi- 
nates the secondary emission, is not sensitive to fluctua- 
tions in the accretion rate. 

In § [5] we examine the validity of the relation f3 corc ~ 
?7p/100. We determine the value of /3 core /%,, its scat- 
ter and its dependence on cluster mass using a sim- 
ple toy model for the evolution of cluster CRs. In or- 
der to describe the effects of mergers, we construct (in 
§ 15.11) the merger hist ory of clusters using the scheme 
of lLacev fc Cold Jl993) , and construct (in § I5.2[) a sim- 
ple model describing the shocks and the adiabatic com- 
pression/expansion of the gas induced by mergers, which 
captures the main features of 3D nu merical merger sim- 
ulations (e.g. iMcCarthv et alJ 12007). Using the simple 
model of § 15.11 and § 15.21 for the evolution of ICM gas, 
we construct in § 15.31 a model for the evolution of the 
CR population. We assume that the effects of diffusion 
are small, and that CRs are advected with the gas. The 
evolution of the CR population is followed taking into 
account energy losses due to Coulomb and inelastic nu- 
clear collisions and assuming that the relativistic parti- 
cles beh aye as an ideal gas w ith an adiabatic index of 4/3 
(see e.g. lEnfilin et al"1l2007l ). The results of this model 
slightly modify out crude estimation to /3 corc ~ ?7p/200, 
with small scatter and weak dependence on cluster mass. 

One limitation of the toy model for ICM evolution is 
that it includes a description of entropy changes that are 
driven only by gravity (accretion and merger shocks). It 



does not include a description of a possible increase of 
the entropy at high redshift by non-gravitational pro- 
cesses (such as s upernovae, st ar formation and galac- 
tic winds, e.g. iKaiserl Il991t lEvrard fc Henrvi 119911: 
iNavarro et all 119951: ICavaliere et alJl!997t iBalogh et al.1 
Il999t iPonman et al.1 Il999f ) , and of the decrease of en- 
tro py due to cool i ng at the centers of ma ssive clusters (see 
e.g. lFabiarjri994 IVoigt fc Fabianll2004l) . As explained in 
§ 15.41 for massive clusters, M > 1O 14,5 M0, these effects 
are unlikely to modify significantly the main conclusion 
of § 15.31 /3 C orc — f? P /200, and are therefore unlikely to 
significantly modify the predicted properties of the non- 
thermal emission. 

In § [5] we explain how our model predictions may be 
tested, and how the values of the parameters rj p and 
r] e may be calibrated, using a controlled sample of clus- 
ters observed in radio and HXRs. We also discuss our 
results' implications for future 7-ray observations with 
space-borne and ground based telescopes. In particular, 
we identify the clusters which are the best candidates for 
detection in 7-rays, and the best candidates for detection 
of emission from pion decays. 

Throughout, a ACDM cosmological model is assumed 
with a Hubble constant Ho = 7O/170 kms -1 Mpc -1 , 
n m = 0.23, n b = 0.039, fi A = 1 - H m , and a 8 = 0.9 

2. A SIMPLE ANALYTIC MODEL FOR THE 
NONTHERMAL EMISSION 

In this section we derive the spectral and radial distri- 
bution of the nonthermal emission produced by ICM CRs 
in massive clusters, M > 10 14 ' 5 Mq. We assume that the 
fraction /3 COTe of plasma energy carried by CR protons at 
the central regions of clusters, which dominate the emis- 
sion from secondary particles, is nearly independent of 
cluster mass and has a small scatter. We show below 
that the secondary nonthermal emission of the clusters 
is determined by their thermal X-ray emission properties 
and by the value of (3 core (the secondary e^ emission de- 
pends also on the strength of the magnetic field within 
the core). The primary nonthermal emission depends on 
the thermal X-ray properties and on r\ e . This allows us 
to derive simple analytic expressions for the spectral and 
spatial distribution of the nonthermal emission. The an- 
alytic expressions show explicitly the dependence of the 
nonthermal emission on model parameters (e.g., the in- 
fluence of a scatter in /3 CO re is easily inferred) . 

Before going into the details of the model, let us briefly 
explain, using si mple arguments, why w e expect /3 C ore — 
rfy/10 (see also iPfrommer et al.l 120071 : Uubelgas et al.1 
l2008f) . As we show in section § [3J the generation of CRs 
is dominated by accretion shocks. Since accretion shocks 
are characterized by high Mach numbers, we expect them 
to produce a flat CR energy spectrum, 



dncn/de cx e 



fl 



For this distribution, the post-shock energy density of 
CR protons within a logarithmic proton energy interval 
(around e) is related to the fraction rj p of the post-shock 
thermal plasma energy deposited in CR protons by 



_ e 2 dn C R, p /de 



Vp 



= gas 



ln[p max /(lGeV/e)] 



(2) 



Here £ gas is the post-shock thermal plasma energy 
density and p max is the maximal momentum of the 
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CR proton spectrum respectively. The maximal mo- 
mentum is determined by comparing the acceleration 
ti me r L , P c/z 4 ~ 2.9 • IOS^-tTi) -1 yr (where v sh = 
8y3T/2mp/3 is the shock velocity, rL, P = 10 2 77/£?-7pc 
is the Larmor radius of the proton, 77 = 7/IO 7 is the 
Lorentz factor of the proton, B-j = B/O.lpG and 
T\ = T/lOkeV is the cluster temperature) to the typ- 
ical cluster dynamical time, idyn ~ 10 9 yr, which gives 
Pmax ~ 10 18 eV/c. Since ln[p max /(lGeV/c)] ~ 20, 
/3cr,p — ??p/20 right behind the accretion shock. 

As we show in section § [SJ merger shocks do not sig- 
nificantly affect the CR population. This, and the long 
(inelastic nuclear collision) cooling time of CR protons, 

*«> - - 2.6 • 10 10 ( 10 _ 3 n cm ^ )" 1 yr, 0) 

where cp™ 1 — 40 mb and n is the gas number density, im- 
plies that the energy density of CR protons produced by 
the accretion shock is later affected mainly by adiabatic 
expansion and compression. The difference between the 
adiabatic indices of the relativistic CRs (4/3) and non- 

relativistic thermal plasma (5/3) implies /3cr,p oc p s &/ 3 , 
where p gas is the gas density. The ratio between the 
mean gas density in the core and the gas density behind 
the accretion shock is typically ~ 10 2 , implying that the 
value of /?cr p at the cluster core is typically expected to 
be ~ ri P /100. 

We have ignored in the preceding discussion the fact 
that different mass elements were accreted at different 
redshifts, hence at a different densities, which implies 
that their density compression factor is not given directly 
by the ratio of densities at the cluster core and at the ac- 
cretion shock. Moreover, major cluster mergers mix the 
gas and heat it with weak shocks, such that some devi- 
ations from the proposed simple relation are expected. 
We present in § [5] a simple toy model for the evolution 
of cluster CRs, which examines the significance of these 
effects. The model yields /3 CO rc — ^ P /200, nearly inde- 
pendent of cluster mass and with small scatter among 
clusters of given mass. We also discuss in §0 the effects 
of entropy increase at high redshift (by non-gravitational 
processes) and of entropy decrease due to cooling at the 
centers of massive clusters. We argue that these effects 
are unlikely to modify significantly the predicted prop- 
erties of the nonthermal emission from massive clusters. 

We derive below analytic expressions for the luminos- 
ity and surface brightness produced by different emission 
mechanisms. For convenience, the details of the calcu- 
lations are given in § \K\ and only the main results are 
given in § !2.2f -§ 12.51 We note that since the emission 
from hadronic processes is concentrated at the cluster 
core, while the emission from primary electrons origi- 
nates from a thin layer around the accretion shock, the 
surface brightness profiles produced by primary and sec- 
ondary particles are different. One may therefore distin- 
guish between the two sources of nonthermal radiation 
by examining the radial dependence of the nonthermal 
emission surface brightness. We first describe in § 12.11 
the thermal emission properties of galaxy clusters, based 
on which the nonthermal emission properties are later 
derived. In § 12.21 and § 12.31 we summarize our model 
results for the emission from neutral (ir ) and charged 
(e^) products of inelastic nuclear collisions of CR pro- 



tons. In § 12.41 we summarize our model results for the 
(IC and synchrotron) emission fr om primary CR elec - 
trons. Motivated by the results of iKushnir et al.l (2009), 
who cons train /3 COTe based on radio em ission, and by the 
results of IKushnir k Waxmanl ((2010) , who constrain r\ e 
based on HXR observations, numerical results are given 
using a normalization of /? corc = 10~ 4 and r\ e = 0.01. 

The analytic results of § I2.2K 12.41 are obtained by us- 
ing simple approximations for the spectrum of secon- 
daries produced and for the IC and synchrotron spec- 
tra (the details are given in § IA"|) . In § 12.51 we compare 
the results of § I2.2I -S 12.41 to those obtained using a more 
accurate parame trization of the spe ctrum of p-p secon- 
daries ( following iKamae et al.l 12006). and using the ex- 
act formulae for IC and sy nchrotron emission (following 
iBlumenthal k Gouldlll970l ). We find that the deviations 
are small. 

The charged secondary emission depends on the the 
strength of the magnetic field in the cluster core. It is 
given in § I2.3l in terms of B/Bqmb, the ratio of the mag- 
netic field to -Bcmb, defined as the magnetic field for 
which the magnetic energy density equals the CMB en- 
ergy density, 

Bcmb = {8naTS MB ) 1/2 = 3.2(1 + zf /iG. (4) 

Since the secondary emission is dominated at high en- 
ergy, > 1 eV, by secondary tt° decay and by IC emission 
of primary electrons, the uncertainty in the value of B 
affects only the predicted radio emission, which is not the 
main focus of the current paper. The deta iled discussion 
of rad io emission from clusters given in IKushnir et al.l 
(2009) indicates that B is within the range of ~ 1/iG to 
~ 10/xG. 

2.1. Clusters of galaxies: thermal emission properties 

The thermal X-ray emission of clusters is usually char- 
acterized by 4 parameters: Lx, the bolometric lumi- 
nosity integrated out to some radius rx from the clus- 
ter center; T, the cluster temperature (weighted aver- 
aged over the cluster emission); f3 and r c , which de- 
fine the radial surface brightness profile, and hence the 
radial density profile j Cavalicrc & Fusco-Fcmianol Il976l : 
iGorenstein et alj|1978t Pones k Formanlll984l ). 

/ r 2X"(3/2),3 
PgasM = pQ ( 1 + -3- J . (5) 

Using the HIFL UCGCS sample of the X-ray -brightest 
galaxy clusters of iReiprich k Bohringerl ()2002l ) , which is 
based on the ROSAT All-Sky X-Ray Survey, we derive 
the following correlations: 

Lx = Lxoh 7 QT" L , 
r c = r- c0 /i 70 1 T 1 Qr , (6) 

where L X q = 2.76-10 45 ergs" 1 , a L = 2.56, r c0 = 223 kpc 
and a r — 1.32. (Note, that various a uthors obtain val- 
ues of (Xl in the range ~ 2.5 — 3, e.g. lM"arkevitchlll998t 
lArnaud k Evrardlll999t IReiprich k Bohringerll2002D . 

In what follows, we give expressions for the nonthermal 
emission as function of T, f3, r c , and Lx- We also give 
expressions which depend on T and /3 only, using the cor- 
relations given in eq. The following cautionary note 
should be made in this context. Our demonstration in 
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§[5l that /3 corc ~ ?7p/200 with weak dependence on cluster 
mass, is based on a model that does not include entropy 
changes not induced by LSS shocks. As explained there, 
and also in § [U early (high redshift) entropy injection 
(which does not significantly affect massive clusters) may 
modify the value of /3 C ore for low mass clusters, thus in- 
troducing at low T an additional, implicit, dependence 
on T through fi cmc . Note, that such early entropy in- 
jections are believed to be responsible for the deviations 
at low T of the observed correlations, eq. ([S]), from the 
self-similar relations expected when entropy changes not 
induced by LSS shocks are neglected (e.g. Lx oc T 2 , see 
lArnaud fc Evrardlll999f ). 

We define the cluster mass, M 2 oo, as the mass con- 
tained within a radius r200; within which the mean den- 
sity is 200 times the critical density, p cl -it. This mass may 
be related to T and j3 by assuming the ICM plasma to be 
isothermal and in hydrostatic equilibrium. Under these 
assumptions, the gravitational mass is given by 



M(< r) 



3/3T 



p,m p G r 2 



(7) 



(the mean molecular weight is fj, ~ 0.59 for fully ionized 
plasma with hydrogen mass fraction x = 0.75), and we 
obtain (assuming r2oo 3> r c) 

/800tt \" 1/2 ( 3/3T V /2 
, 200 ^__ Pcrit j 

= 3.1/3 1/2 T 1 1/2 /if 1 Mpc, 

/800tt V 1/2 / 3/3T \ 3/2 

= 3.5-10 15 ^ 3/2 T 1 3/2 /i 70 1 Af Q . (8) 

In addition to the assumption that the ICM is isother- 
mal and in hydrostatic equilibrium, we have extrapolated 
here the surface brightness given by eq. ([5]) out to r2oo ? 
which may be larger than rx- This implies p oc r -2 at 
large radii, in contras t with the p <x r~ 3 d ependence ex- 
pected at large r fe. g.lNavarro et aLlll997fl . The d etailed 
discussion given m Rci prich fc Bohringerl (|200l of the 
accuracy of cluster mass determination under these ap- 
proximations shows that eq. (jSJ may overestimate M200 
by no more than 20%. 



2.2. 



decay emission 



In this section we present our model results for the 
emission from 7r° decays. Since the CR proton spec- 
trum is flat, the p-p 7-ray luminosity per logarithmic 
photon energy bin, vLf, p , is energy independent above 
~ 0.1 GeV (the threshold energy for pion production 
is e t h — 1.22 GeV and the photon energy is ~ 0.1 
of the parent CR proton energy). The observed spec- 
trum is expected to be cut off at high energy due to 
pair production interactions of the high energy photons 
with infra-red background photons. For nearby clus- 
ters, the 77 — > e+e~ c utoff is expected at ~ lOTeV 
(|Franceschini et aT1l2008|) . The tt° luminosity is given by 
(see eq. (|Al3]) ) 



i/LPP~1.3-10 41 ) 8 corei _4r 1 J 



1/2 



L 



x 



h~ 2 3 ■ 10 45 ergs" 1 



ergs 



(9) 



where /3 C ore,-4 = /3 C ore/10~ 4 . vL'PJ' depends lin- 
early on Lx, since both hadronic emission and ther- 
mal bremsstrahlung emi ssion depend on density squared 
(|Katz fc Waxmar] [2008). Using the correlations eq. ([6]) 
we have 

i/LgP ~ 1.2 • lO^ftore^Tf 06 ergs- 1 . (10) 

The p-p 7-ray surface brightness above some energy 
£iy,min at some distance f = r /r c f rom the cluster center is 
given by (for /3 > 0.5, see eq. (|A5|) and lSarazin fc Bahcalll 
I1977T) 



^ min (f)^5.5-10- 7 Up -0/3, 



core,— 4^1 



1/2 



■>x 



/if 2 3 • 10 45 erg s- 1 / V ^ro 200 k P c 

max(£^ min ,0.l£ t h) 
10 GeV 

f-i , -2\ -3/3 + 1/2 _ 2 -1 -1 

x(l + r^) phem s sr x . 

Using the correlations eq. ([6]) we have 

SJP (f)~4.1 • 10- 7 Up - I ) /3 corc ,_ 4 T 1 1 



(11) 



max e. 



2, 
O.leth) 



x (1 + f* 



10 GeV 

n -3,3+1/2 



phcm-^s-^i'- 1 . (12) 

Note, that in the derivation of the surface brightness 
given above we have neglected the variation of /3cr,p with 
radius, and used /3cr,p = /?coro- This is justified since 

the surface brightness S is oc /3cR,pPgaf 1 ^ 3 ' 3 , and, as 

1/3 

we show in § [5j /3cr,p oc p gas ■ This implies that the 
variation of the surface brightness is dominated by the 

Pgaf 1 ^ 3 ' 9 term (for typical /3 values, which are in the 
range of 0.5 to 1). 

2.3. Emission from charged secondaries 

p-p collisions also produce secondary electrons and 
positrons, which cool by emitting synchrotron radiation 
and by IC scattering of CMB photons. Since variations 
in the secondary injection and in the cluster magnetic 
field occur on time scales similar to, or larger than, 
the dynamical time scale of the cluster, idyn ~ 1 Gyr, 
the distribution of secondaries is in a steady state at 
high energies, at which the secondaries lose all their en- 
ergy to radiation on a time scale short compared to the 
cluster dynamical time. For magnetic field values of 
1 — 10 /iG, the Lorentz factor of secondaries with cooling 
time, Qirm e c/ {B 2 + -B 2 , mb )ot7, which equals the dynam- 
ical time, is 200 < 7 cool < 2000. 

The luminosity due to IC scattering of CMB photons 
by secondaries is given by (see eq. (|A8|I ) 



uL]p' e -3.3-10 



40 



d2 

-"CMB 



R2 

CMB 

Lx 



B 2 



core,— 4-^1 



h m z 3 ■ 10 45 ergs- 1 



ergs 



1/2 



(13) 



Eq. (fT3"|) holds for photon energies in the range 
7c 2 oo1 3T cmb(1 + zf < e ph < 7 2 iax 3T CM B(l + zf, where 



6 



7max — 0.1e max /m e c 2 (e max is the maximal energy of the 
CR protons, and the secondary energy is ~ 0.1 of the par- 
ent proton energy). The surface brightness is given by 
(see eq. (|A"9t ) 



X /3corc,— 4?\ ^ ( 



R2 

D CMB 



# 2 4- R 2 
CMB ' D 



L 



x 



/if 2 3 • 10 45 ergs- 1 
2 _ 1 

10 GeV/ 



ftfo 200 kpc 

x (1 + f 2 )~ 3/m/2 phcm- 2 s- 1 sr- 1 . (14) 
Using the correlations eq. ((6|) we have 

vL IC > e± ~<i n • 1 o 40 B ™ B 

iU R2 i R2 

"CMB ' -° 



and 



x ^ coro ,_ 4 T 1 3 - 06 ergs" 1 
5^1 (f)~ 1.0 -10- 7 (3/5-0 



(15) 



R2 

C CMB 



^CMB + ^ 2 



v o rpO.42 / £ ^,min \ 1 

x (1 + r 2 )" 3/?+1/2 phcm- 2 s- 1 sr- 1 . (16) 
Note that although the 7-ray luminosity produced by 
7T° decays is bigger above ~ 0.1 GeV by a factor 
~ 4(1 + B 2 /Bq MB ) than that produced by charged sec- 
ondaries (compare eq. © and eq. (1731) ). the secondaries' 
emission is the main hadronic emission mechanism below 
this energy. 

For photons in the energy range 7 2 00 i£o < £ph < 
7max £ o, where £0 = ZfaeB /2m e c, the seco ndary syn- 
chrotron luminosity from is given by (see eq. (|A12|) ) 

V LT^ e± ^3.3 ■ 10 4 ° „o B2 ^ flcorc,-^ 2 



r2 1 R2 ' 

CMB "T" 

T ->x 



I erg s 

/if 2 3- 10 45 ergs-V 6 ' 



(17) 



and the surface brightness is (see eq. (|A13jl ) 



X /3corc,— 4-^1 



1/2 



A" 



h™ 1 200 kpc 



/if 2 3- 10 45 ergs" 1 
V 1.4 GHz/ 



x (1 + f 2 ) 3/?+1/2 mJy arcmin" 2 . (18) 



Using the correlations eq. (JB]) we have 

B 2 



40 



^CMB + B 2 



xT 3.06 ergs -l 5 



(19) 



and 



S r c,e± (f) ^ 9 . 8 .( 3/3 _3^__& 



9 / R 2 4- R2 
z / CMB ' 



x /3 coro _ 4 T° 42 f 

H ' 1 V 1.4 GHz/ 

x (1 + f 2 y 3P+1/2 mJy arcmm -2 . (20) 



Note, that due to the uncertainties in B and t^ yn , 
which lead to uncertainty in 7 CO oi, our estimates of the 
emission from secondary electrons and positrons with 
7 ~ 7cooi are uncertain. 

2.4. Primary electron emission 

We now turn to estimating the radiation produced by 
primary electrons at the accretion shock. We assume 
that the clusters are spherical and that gas is accreted 
onto clusters in spherically symmetric flows at a rate 

AI = /inst-^ooA-F/- /inst is a dimensionless parame- 
ter of order unity, reflecting the temporal fluctuations 
of M j (M2oo/ij?)- As discussed in §[3l 3D numerical sim- 
ulations indicate that the average value of /i ns t is rs 0.5. 

In order to calculate the surface brightness of the pri- 
mary electrons' radiation, it is necessary to estimate 
the position of the accretion shock. We assume here 
that the accretion shock is located at r ~ ^200 j since 
spherical collapse models predict a cluster virial den- 
sity < p v i r >~ 178p cr it (for fi m = 1 and J7a = 0, 
with weak dependance on the background cosmology for 
0.3 < f2 OT < 1). The validity of the simplifying assump- 
tions described above is tested in §[3]by comparing our re- 
sults to those of numerical 3D simulations. We find that 
the results obtained in this section are in good agreement 
with the results of detailed numerical simulations. 

As in the case of the charged secondaries, we assume 
that the distribution of the primaries is in a steady 
state, since electrons (at the relevant energies) lose all 
their energy to radiation on a time scale short com- 
pared to the cluster dynamical time, i^yn ~ 1 Gyr. Un- 
like the secondary electrons and positrons, which lose 
energy through both IC and synchrotron emission, pri- 
mary electrons lose their energy mainly by IC scattering 
of CMB photons, since the magnetic field at the accre- 
tion shock is expected to be weak, ~ 0.1 fiG <C -Bcmb 
(|Waxman fc Loebl feOOO) . Thus the steady state assump- 
tion holds for primaries with Lorenz factors 7 > 7 CO oi ~ 
2000. The maximal energy of the primary electrons is de- 
termined by equating the acceleration time, r^^ c c/v 2 h ~ 
1.6-10 3 77(-B-7Ti) _1 yr (where r L , c = 5.5-10" 2 7 7 /S_ 7 pc 
is the Larmor radius of the electron), to the IC cooling 
time, 67rm e c/i3Q MB o'T7 — 2.3-10 5 7^~ (l + z)~ 4 yr, which 
yields 7max ~ 1.2 • 10 8 /bIt7Y(1 + z)- 2 . 

For photon energies 7 2 oo1 3T C mb(1 + z) 4 < e ph < 
7max3TcMB(l + z) 4 , the primary I C lum inosity is give n, 
to very good accuracy, by (see eq. (|A16I) and eq. (|A17D ) 

^IC, sh ock^ 1 8 . 10 43 ( /inst77e )_ 2/3 3/2 

A-^T* /2 Z(z)erg S -\ (21) 

Here, (/mst7? e )_ 2 = /mstWl0~ 2 , fb = ^b/^m and 
Z(z) = (tHH(z))^ 1 (Z(z) depends weakly on the as- 
sumed cosmology). Assuming the emission originates 
from a thin layer with thickness w behind the accretion 
shock, the surface brightness is given by 

S»f(r) ~ 1.5 • 10- 7 (f inst r) e )_ 2 fi 1 ' 2 



fb 



11 VToc 



Q.17)~ L V 10 GeV/ 
x £(r/r 2 oo,w/ r 20o) 
x Z(z)h 2 (z) phcm - ' 4 s 



X 1 



- 2 -- l Bt-\ 



(22) 
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where 



3(yi^-yJ(l-y) 2 -x^ 



i-(i-a) J 



x ^ 1 -y (23) 

1 - y < x < 1. 



For w <C ?°2oo, ^(x, y) can be approximated in the regime 
x < 1 - y by 



1 



vT 



(24) 



The thickness w of the emitting region is approximately 
given by the product of the cooling time of the emitting 
electrons, i C ooi, and the velocity of the downstream fluid 
relative to the shock velocity, Ud- Since idyn ~ ^oo/wd, 
the approximation w <C ^200 holds for 7 > 7 coo i ~ 2000. 
Note that the rapid increase of the surface brightness 
near the accretion shock, inferred from eq. (|24|) . is likely 
to be suppressed by small deviations from spherical sym- 
metry. 

Below ~ 1 TeV, the primary IC 7-ray luminos- 
ity is larger than the secondary 7-ray luminosity by 

a factor ~ 150(/ ins t77 e )^ 2 (^corc)lI/3 3/2 r i _1/2 (compare 
eq. (|2T|) to the neutral component contribution eq. (|T0|) . 
which dominates the secondary emission at this en- 
ergy range). The primary IC HXR luminosity is larger 
than the secondary HXR luminosity by a factor ~ 

600(/ inst r ?e )_ 2 (/3 corc ):i/3 3 / 2 T 1 " 1/2 (l + B 2 /B 2 CMB ) (com- 
pare eq. (I21[) to the charged component contribution 
eq. (|15p , which dominates the secondary emission at this 
energy range). However, since the primary electrons' IC 
surface brightness is lowest at the cluster core, while the 
secondaries' emission surface brightness follows the ther- 
mal surface brightness (see e.g. eq. (|TT|) ). the surface 
brightness due to these two sources of radiation may be 
comparable at the cluster core. 

For photon energies 7 2 ool £ < £ P h < 7ma X £ o, the pri- 
mary electrons' synchrotron luminosity is given by (see 
eq. jXM ) 

^r C,Sh ° Ck - 1-7 • 10 40 (finstVe)^ /3 3/2 
fb \ rp5/2j^2 



0.17 y 

x Z(z)(l + z)~ 4 ergs -1 . 



(25) 



This luminosity is comparable to the luminosity pro- 
duced by secondaries. However, the primary electron 
synchrotron surface brightness, 

5 syac,shock (r) „ ^ 4 . 1Q -2 (f^J^ 
fb 



T 3/2 / 

11 Via 



s 0.17y 1 V 1.4 GHz, 

x B -ii (r/r 2 oo,w/r 2 oo) Z(z) 

x ft 2 (z)(l + z)~ 4 mjyarcmin~ 2 (26) 

(see eq. JSHJ)), is negligible compared to that produced 
by secondaries (see eq. ([15])). 

2.5. Summary of results 

Figure Q] shows the nonthermal luminosity, vL vl as 
function of photon energy for a T = lOkeV, /3 = 2/3 
cluster (with Lx and r c determined by eq. ©), with 
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Fig. 1. — vh v as function of photon energy for a T = lOkeV, 
= 2/3 cluster, with p coTC = 10~ 4 , fmstVe = 0.01, B = Bomb at 
the cluster core and B = 0.1 /iG at the accretion shock. Solid lines 
show the simple estimates given by eqs. 10, Ill3t . 117t . (1 2 1 1 and 
I l25t . a nd dashed li nes are obtained using the parametrization of 
IKamae ct al. ( 2006|) for the secondary spectrum of p-p interactions, 
and the exact formulae for IC scatter ing and synchrotron emission 
given by Blumcnthal & Gould (1970). Suppression of the flux by 
pair production (in interactions with IR background photons) was 
included, applying a suppression f actor exp(— r-yy) with pa ir pro- 
duction optical depth (taken from Franccschini ct al. 2008) corre- 
sponding to a cluster at the distance of Coma. Note, that there 
is an artificial high energy cut-off in the sec ondary emission spec - 
tra calculated using the parametrization of Kamac et al.l (2006), 
resulting from the upper limit on the proton energy used in this 
parametrization, 512 TeV. 



j3 corc = 10~ 4 , /inst^e = 0.01, B = B C MB at the cluster 
core and B = 0.1 fj,G at the accretion shock (for our as- 
sumed background cosmology we have Z(z = 0) ~ 0.97). 
Solid lines show the simple estimates given above by 
eqs. ©, dHJ), QUI, (@TJ and ([25]), and th e dashed lines 
are ob tained using the parametrization of Kam ae~et al.l 
(2006) for the secondary spectrum of p-p interactions, 
and the exact for mulae for IC scattering and s ynchrotron 
emission given bv lBlumenthal fc Gouldl (fl970h ■ Suppres- 
sion of the flux by pair production (in interactions with 
IR background photons) was included, applying a sup- 
pression factor ex p(— t^) with pair produ ction optical 
depth (taken from iFranceschini et al.ll2008f ) correspond- 
ing to a cluster at the distance of Coma. Note, that 
there is an artificial high energy cut-off in the secondary 
emission spec t ra cal culated using the parametrization of 
IKamae et al.1 (|2006l ). resulting from the upper limit on 
the proton energy used in this parametrization, 512 TeV. 
However, the secondaries' emission is always small in this 
range compared to other emission processes. 

In fig. [2] we show HXR, 7-ray and radio surface bright- 
ness profiles for the same cluster at various energies 
(see eqs. (HJ), QU), (18]), © and (|26])V The thick- 
ness of the layer behind the accretion shock, from which 
the primary electron emission originates, is estimated 
to be w ~ Udtcooi- For the chosen cluster parame- 
ters, the dominant nonthermal emission process in the 
HXR (> 100 keV) band is IC emission from primary elec- 
trons at the accretion shock. At higher photon energies 
(> lOGeV) the emission from pion decays at the core 
of the cluster becomes comparable to the IC emission 
from the accretion shock electrons. Note, that the IC 
emission from accretion shock electrons is uncertain at 



S (>100keV) S (>10GeV) 




S (>1 TeV) S (@ 1.4GHz) 




Fig. 2. — HXR, 7-ray and r adio s urfa ce br ightne ss a s function 
of radius, given by eqs. JTTJ, ITU ). (Tl8l . ll22l and (f26l . Cluster 
parameters and line types are the same as in fig. [l] The thickness 
of the layer behind the accretion shock, from which the primary 
electron emission originates, is estimated to be w ~ u,jt coo i <S r 200- 
Note that the rapid increase of the surface brightness near the 
accretion shock is likely to be suppressed by small deviations from 
spherical symmetry. 



energies exceeding ~ 1 TeV, due to the uncertainty in 
the maximal energy, ~ few TeV, to which electrons may 
be accelerated in such shocks. Thus, although at the 
figures shown the IC emission dominates the luminosity 
also at the highest energies, > 1 TeV, pion decays may 
become the dominant emission process at these energies. 
At all energies, the contribution from secondary IC is 
small. However, the synchrotron emission of secondaries 
is the dominant emission process at the radio band. 

3. COMPARISON TO DETAILED NUMERICAL 
SIMULATIONS 

In this section we compare our model's results with 
a variety of detailed numerical simulations. We show 
that our results are in agreement with those of numerical 
simulations, with the exception of deviations which are 
d ue to inaccuracies of the numerical calculations. 

IKeshet et ail (|2003t) used TreeSPH simulations of LSS 
formation to estimate the IC emission from electrons ac- 
celerated in LSS shocks. The typical > 10 GeV f lux 
from rich clusters obtained by Kc shet et al.l (|2003l) is 
few x 10~ 7 (?7 e /0.05) phcm~ 2 s _1 sr _1 (see their fig. 10) 
consis tent with the prediction of eq. ([2"2")l . Keshe t et al.l 
(2003) also used the simulations to determine N(> /), 
the number of sources expected with flux exceeding / (at 
a given photon energy threshold). The model presented 
here for the nonthermal emission of clusters may be used 
to predict iV(> /) for a given number density of halos 
as function of redshift and hal o mass. Such an exercis e 
has already been carried out bv lWaxman fc Loebl (2000), 
who derived analytic expressions for the 7-ray luminosity 
of primary electrons and u sed them to predict N (> /) as- 
suming a Press-Schechter ([Press fc SchechterlH97l ' mass 
function for the cluster halo density. The expressions de- 
rived here for the gamma-ray luminosity of a cluster of 
a given mass a re identical, up to a no rmalization fac- 
tor, to those of IWaxman fc Loeb (2000), implying that 
our prediction fo r N(> f) would be similar to that of 
IWaxman fc Loebl (2000) . ' Since in the normalization used 



bv lWaxman fc Loebl (|2000[ ) the accretion rate is ~ 3/^g t 
larger than the one used here, and since the IC accre- 
tion luminosity is proportional to the accretion rate, the 
luminosity function, N(> /), obtained using our normal- 
ization is related to the one derived bv lWaxman fc Loebl 
S2mh. AW> f) by N(> f) = N WL (> 3/^/). 
IKeshet et al.l (|2003l ) found that their numerical source 
number cou nts, jV n um(> f), fall sho rt of the analytical 
prediction of lWaxman fc Loebl ([2000) by a factor of ~ 6, 
in the sense that N num (> f) w %i(> 6/). This im- 
plies that the analysis presente d here is consi s tent w ith 
the numerical number counts of Kes het et al.l ()2003l ) for 
/inst — 0.5. 

The results of TreeSPH simulations of LSS evolution 
incorporating CR generation in L SS shocks were recently 
reported in a series of papers (|Pfrommer et al.l I2007L 
120081 iPfrommerl [2008h . In these simulations, rj p = 0.5 
and rj e = 0.05 were adopted, and the CR evolution was 
followed under the assumption that the CR energy spec- 
trum is a power-law. Figure [3] presen ts a comparison 
of our results with those obtained by iPfrommer et all 
(120081) for a massive (M = 10 15 Mq, T = 9.6 keV) 
merging cluster, for which detailed emission spectra are 
given. The simulation's surface brightness profiles of 
pion decay and secondary IC emission above 100 MeV 
follow the thermal surface brightness with central values 

of 5P p (0) ~ 5 ■ lO^phcm^s-isr" 1 and S£ c ' e± (0) ~ 
5 ■ 10~ 5 phcm~ 2 s -1 sr -1 respectively. The profiles and 
the central values of the surface brightness of pion decay 
and secondary IC emission obtained in the simulation 
are consistent (to within a factor of two) with those pre- 
dicted by eqs. JT2J and ^TBJ) for /3 core = 2 • 10 ~ 3 and 
P = 2/3. This value of /3 C ore is consistent with our sim- 
ple model prediction, fi cote ~ 7/ p /200, for the value of r/ p 
chosen in the simulations, rj p — 0.5. The ratio of pion 
decay to secondary IC luminosity obtained in the simu- 
lation is larger than predicted by our model. This may 
be a result of using power-law approximations for the 
CR emission spectra (compare the solid and dashed lines 
of our model near 100 MeV in figs, [TJ [3]). We note that 
the numerical calculations, to which our model predic- 
tions are compared in fig. [3l included radiative cooling 
of the ICM plasma, which is not included in our model. 
As explained in § 15.41 such radiative cooling is not ex- 
pected to significantly affect the nonthermal emission of 
massive clusters. This conclusion is consistent with the 
results of the numeric al calculations (see, e.g., fig. 13 in 
IPfrommer et al.ll2008D . 

The value of /3 C orc obtained in the simulation for the 
cluster used for the comparison of fig. [3] can not be eas- 
ily extracted from the results reported, since a profile 
of the ratio of CR to thermal gas pressure is not given. 
Instead, a profile of CR to thermal gas pressure ratio 
averaged over 9 clusters, including the one used for the 
comparison of fig. [3]and 8 smaller (5 • 10 13 — 5 • 10 14 M Q ) 
clusters, is given. This averaged pressure ratio profile 
is nearly constant throughout the cluster at a value of 
~ 0.1, and rises sharply within the inner 2% of the virial 
radius to ~ 1. A pressure ratio of ~ 0.1 corresponds 
(for r] p = 0.5) to /3 CO re/f?p ~ 2/100, larger than inferred 
by our analysis. However, this profile is probably bi- 
ased towards larger values of j3 core /rj p due to the large 
number of small clusters for which radiative cooling is 
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Fig. 3.— Figure 12 of IPfrommer et alj J200|) (curtesy of C. 
Pfrommer). Over plotted are the su rfac e brig htness cur ves ob- 
tained from our simple analysis (eqs. 1111 , AMU and 1221 '), for a 
T = 9.4keV, /3 = 2/3 cluster with /3 covc = 2 ■ 10 -3 , rj e = 0.05 (as 
chosen in the simulation), and /i ns t = 1 (Lx and r c were deter- 
mined according to eq. (6}). /Score = 2 ■ 10 — 3 is consistent with our 
model prediction, /3 C ore ~ ^p/200, for the value of r\ v chosen in the 
simulation, r] p = 0.5. Line types are the same as in fig. [T] Note 
that the rapid increase of the (primary) surface brightness near the 
accretion shock is likely to be suppressed by small deviations from 
spherical symmetry. 



more important (radiative cooling also leads to the sharp 
rise at the inner 2% of the virial radius). By examining 
the pressure ratio profiles of M ~ 10 14 M Q clusters in 
the non-radiative simulations, which should be very sim- 
ila r to those of more ma ssive clusters (see also table 3 
in IPfrommer et al"1l2007l ). we find that indeed this ra- 
tio decreases toward the center of the cluster to a value 
/W?7p ~ 1/500. 

The primary electrons' IC surface brightness obtained 
in the simulation is approximately uniform across the 
cluster. Its value, Sl C ' shock ~ 5 • 1CT 6 phem" 2 s" 1 sr" 1 , 
is lower by a factor of ~ 10 compared to the prediction 
of eq. ([22]) (for ,5 = 2/3 and /i nst = 1)- This discrepancy 
is due to the fact that the accretion shock is located 
at the simulation at a radius of ~ 10 Mpc, ~ 3 times 
larger than r 2 oo (see eq. ©), where we assume the shock 
to be located. That is, the total primary IC luminos- 
ity obtained in the simulation is in agreement with our 
prediction, eq. (|2"T|) , but the surface brightness is ~ 10 
times lower since the shock is located in the simulation 
at a ~ 3 times larger radius. We believe that the shock 
radius is overestimated in the simulation, since our re- 
sults are consi stent with those of the n umerica l simu- 
lations of both iKeshet et ail (|2003l) and iMiniatil (|2003l) 
(see below), which imply the shock position to b e closer 
to r2 no than obtain e d by Pfrommer et al.1 (|2007| ). More- 
over, [Molnar|etIal] (|2009 ) have recently shown that ac- 
cretion shocks are located in SPH simulations at a radius 
which is ~ 3 larger than that obtained in AMR simula- 
tions, probably due to numerical inaccuracies in the SPH 
simulations. 

Figure |H pre sents a g r aphic al comparison of our results 
with those of IMiniatil (|2003f ) . who modelled the evolu- 
tion of CR protons and electrons using an Eulerian+N- 




FlG. 4. — Figure 2 from IMiniatil Q2003T ) (curtesy of F. Miniati). 
Over plott ed a re vL v cu rves obtained from our simple analysis 
(eqs. ©, JTS} and (21)), for a T = 4keV and /3 = 2/3 cluster 
with /Score = 0.05, T) e = 0.01 (as chosen in the simulation), and 
/inst = 1 (Lx and r c were determined according to eq. (|6]|). The 
large of value of /3 C ore obtained in the simulations is probably due 
to the low resolution of the simulation (see text). Line types are 
the same as in fig. [l] 



body code to describe the evolution of the LSS, assuming 
CRs are generated at strong shocks with rj p ~ 0.6 and 
r) e = 0.01 . The IC lum inosity of primary electrons ob- 



tained bv IMiniatil (|2003l ). vL]?> 



hock 



2 • 10 42 ergs~\ is 



consistent with the prediction of eq. (f2"Tj) (for (3 = 2/3 
and /i ns t = 1)- The value of /3 C ore obtained in the simula- 
tion, P corc ~ 0.05 for a T = 4 keV cluster, implies a ratio 
flcorc/Vp ~ 1/10, significantly larger than the value we 
estim ated, ~ 1/200, which is consistent with the results 
of the IPfrommer et al.l (|2007f ) simulations. This discrep- 
ancy is most likely d ue to the low, ~ 100 kpc, resolution 
of the simulation of IMiniatil ((2003) , which doesn't al- 
low one to properly resolve the adiabatic compression 
at the core. The pion decay lumi nosity and the sec- 
ondaries' IC luminosity obtained by IMiniatil ()2003l ) are 

zaLpp ~ 2 • lO^ergs" 1 and vL l ^ e± ~ 5 • lO^ergs" 1 
respectively. Both are comparable to our model predic- 
tions, eqs. (HO} and (TTS]), for T = 4keV and /3 coro = 0.05. 

4. COMPARISON TO OBSERVATIONS - THE COMA 
CLUSTER 

In this section we compare the results of our ana- 
lytic model to various observations of the Coma clus- 
ter. We use these observations to test our model 
and to constrain its parameters. For the compari- 
son, we use the following parameters to characterize 
the Coma cluster: Iboma = 8.25 keV, Lx,Coma = 1-1 • 
10 45 /if 2 ergs~\ /3coma = 0.654, r CjC oma = 246/ifokpc, 
t-200, Coma - 2.3ft-7 1 Mpc, M 2 oo,co ma ^ 1.4 ■ 1Q 15 Mq and 
z Coma = 0.0232 (|Reiprich fc Bohringerl I2002T) . With 
these parameters we have r c ,coroa/<^Coma — 0.14° and 
r200,Coma/rfcoma — 1-3°, where d C oma is the distance 
to the Coma cluste r. Note that the surface b right- 
ness measurements of lReiprich fc Bohringerl (|2002l ) reach 
rx,Coma — 2.9/ifo 1 Mpc > r 2 oo,Coma, so extrapolation is 
not needed in order to determine r2oo and Af2oo for this 



10 



cluster. 

In order to make the comparison with observations 
more straightforward, we first give below explicit expres- 
sions for the flux predicted by our model within a disk of 
angular radius 8 (centered at the cluster's center) for a 
cluster at a distance d ~ cz/H . We further assume that 
within the cluster core B 2 ^> Bn„ AVt , as inferred from 
radio observations (|Kushnir et al.ll2009h . and normalize 
our results to £>_ 5 = 5/10 ^G. Using the results of §0 
we have 



F»>» mi M = 1-4 ' 10" 11 3/3 - - /3 coro ,_ 4 T 1 



-,1/2 



X 



h^3 ■ lO^ergs" 1 



finin(ed/r c ,r 2 oo/i , c) 

X I ( 1 +^ 2 )" 







max 



i,0.1e t h) 



lOGeV 
x h% phcm~ 2 s^ 1 
for 7-rays from pion decay, 

1-13 



-l 



■3/3+1/2 



rdr 

-2 



(27) 



3.4-10 ::! f 3/3 — | ) B r (,.,,,... , 



-■1/2 



Lx 



x T 

V/iTO 3 - 1045 ergs" 1 



'70 ' 

min(6>d/r c ,r2oo/>"c) 



(l + r 2 ) 3 ^ +1 / 2 r( ^ r 



VIOGeV/ 



x hj ph cm 2 s 1 
for IC emission from secondaries, 

S s J nc - e± (8) = 3.7 ^3/3 - /3 coro ,-4T 1 1/2 

Lx 



(28) 



/if 2 3- 10 45 ergs- 1 

/■min(ed/r c ,r2oo/' , c) 

< / ( 1 + f2 ) 




2 \ -3/3+1/2 



~'dr 



x /i2 Jy 



.1.4 GHz/ V-'-cmu.-i. 

l 70°y (29) 
for synchrotron emission from secondaries, and 

^>Sr; k w =4.7 • io- 7 «u^)_ a /3 1 / 2 

/ fb \ ^3/2 / giynin \ _1 

*Ui7/ 1 VlOGevJ 

x g acc .(8)Z(z)h 2 w (z) phcm- 2 s- 1 , (30) 



with 



ffacc.(^) = 2#2 



1" , (31) 

"200, 



for IC emission of primary electrons. Here #200 — ^200/^ 
and (f hist) 9 is the average value of /i ns t over the disk 



considered. Eq (|31[) gives an approximate description of 
the dependence of F on 8, for the case where w -C r2oo 
(see eq. (I2~4"10 . For w = 0.1r 2 oo, eq. (j3"Tj) is accurate to 
better than ~ 25%. 

4.1. HXR observations 

An excess of HXR emission over the expected 
thermal bremsstrahhmg emission has been ob- 
served in the Coma cluster with instruments 
on board three dif f erent X-ray satelli t es: RXTE 
(iRephaeli etail Il999t IRephaeli fc Gruberl I2002D. Bep- 
poSa x |Fusco-Femiano et al.lll999H Fusco-Fcmiano et al.1 
2001 120071) and IN TEGRAL (jEckert et all [20071 
Lutovinov et al.l 120081 ) . Since the observations are in 
agreement with each other, we focus on the recent 
INTEGRAL observations. The INTEGRAL-measured 
flux in the 44 — 107 keV band (where the thermal 
contribution is small) is (1.8 ± 1.1) • 10~ n ergcm~ 2 s -1 
within 8 < 1°. To study the spatial structure of the 
HXR emission, images in the "soft" 17 — 28.5 keV and in 
the "hard" 44 - 107 keV INTEGRAL bands were used. 
In the soft band, Coma is clearly an extended source. 
At the hard-band, the raw image does not show any 
significant substructure or correlation with the cluster's 
thermal emission (on 1° scale). The extended nature of 
the HXR emission implies that the radiating particles 
are not secondary particles produced in the interaction 
of cosmic-rays with the ICM, since the generation of 
such secondary particles should be strongly concentrated 
towards the cluster's center. 

In fig. Owe show the estimated HXR flux within a disk 
of angular radius 8 as function of 8 for the Coma clus- 
ter in the 44 — 107 keV band, assuming /i ns t?7e = 0.01, 
/3 corc = 10" 4 and B = 10 fxG. The INTEGRAL mea- 
surement and the thermal emission in this band are also 
given. Since the nonthermal flux is dominated in this 
band, according to our model, by IC emission of accre- 
tion shock electrons, the observed HXR flux may be used 
to calibrate rj e fi nsi (/3 core is not constrained since the con- 
tribution from secondaries is negligible). According to 
our model, the hard-band image should not show any 
significant substructure since it originates from the ac- 
cretion shock, and no spatial correlation with the cluster 
thermal emission is expected. 

Using eqs. (QHD and flUJ), with Z(0) ~ 0.96 and 



1°/0.14° 



(1 



-2\-3/3com a + l/2 



fdf ~ 0.90, 



(32) 



we have 



-Fcoma,total(l°) — 6.7 ■ 10 16 /3 C orc, -±B_ 2 



erg 



+ 1.4-Hr^((/ inst ) 1 o7 7e )_ 2 — f-(33) 

cm z s 

(where we have multiplied the number flux by 

ln(£max/£min)(l/£min - with £ min = 44 keV 

and £ max = 107 keV, in order to obtain the energy 
flux). The INTEGRAL measurement thus constrains 
0.06 < (/inst)i°f?e < 0.23. Note, that our analysis also 
rules out secondary emission as the source of the HXR 
excess, since in order to reproduce the observed flux by 
secondary emission /3 corc ~ 5 would be required, imply- 
ing that the CR energy density is much greater than the 
thermal energy density of the ICM. 
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Fig. 5. — The predicted HXR flux within a disk of angular radius 
9 as function of 6 for the Coma cluster in the 44 — 107 keV band, 
assuming /i ns t»?e = 0.01, /3 COTe = 10~ 4 and B = lOttG. Line types 
are the same as in fig. [TJ The error-bar represents INTEGRAL'S 
measurement. According to our model, the flux is dominated by 
IC emission of accretion shock electrons, and INTEGRAL'S mea- 
surement implies 0.06 < (/i ns t)l or ?e < 0.23. The HXR excess mea- 
sured by INTEGRAL can not be due to secondary emission, since 
/3 C orc ~ 5 would be required in order to reproduce the observed 
flux by secondary emission, implying that the CR energy density 
is much greater than the thermal energy density of the ICM. 



4.2. EGRET's "/-ray observations 

iReimer et all (|2003| ) report EGRET upper limits on 
high-energy 7-ray emission from the Coma cluster. 
EGRET's upper limit for the flux above 100 MeV is 
Fcoma < 3.81 • 10~ 8 phcm~ 2 s _1 (with Coma assumed 
to be a point source; EGRET's FWHM at this energy 
is 5.8°). In fig. [6] we compare the estimated 7-ray flux 
within a disk of angular radius 9 as function of 9 for 
the Coma cluster with an energy threshold of 100 MeV, 
assuming / ins t?7e = 0.01, /3 cor c = 10~ 4 and B = 10 //G, 
with EGRET's upper limit. Since the flux is dominated 
in this band, according to our model, by IC emission of 
accretion shock electrons, EGRET's upper limit sets an 
upper bound for rj e f msi (and only a weak constraint on 
/3 core ; Note that the tt° decays contribute significantly 
only above ~ 122 MeV due to the threshold for pion pro- 
duction) . 

Using eqs. (27]), (@SJ and flU, and 



u 



(1+f 2 ) 



-3/3 C oma + l/2 



fdf~ 0.95, (34) 



we have 

^Coma, total - 



:7.0- 10" n /3 c 



+ 3.0- 10- 8 ((/ in8 t>i. 3 °»7e)_ 



ph 



(35) 



(since the energy band is close to the pion produc- 
tion threshold, we used the more detailed spectral de- 
pendence of the secondaries, as described in § 12.51 to 
obtain F Coma ,pp — 7.0 • 10 _11 ^ CO ro,-4 phcm~ 2 s _1 in- 
stead of Fcoma.pp ^ 2.0 • 10~ 10 /3 corc _ 4 phcm" 2 s" 1 ). 
EGRET's upper limit implies therefore /3 core < 0.05 and 
(/inst)i.3°?7e < 1.4 ■ 10~ 2 . We don't regard the upper 
limit on fi ns tVe to be in contradiction with the range 
calibrated in § 14. 1[ since it may be explained by a spa- 
tial dependence of /i ns t, or by a primary electron energy 
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Fig. 6. — The predicted > 100 MeV 7-ray flux within a disk of 
angular radius 9 as function of 8 for the Coma cluster, assuming 
/inst^e = 0.01, /3 corc = 10 -4 and B = 10 fiG. Line types are the 
same as in fig.[T] The arrow represents EGRET's upper limit. The 
upper limit implies (/i ns t)l.3° Ve < 1-4 • 10~ 2 and /3 C orc < 0.05. 



distribution slightly steeper than dn/de oc e~ 2 . Future 
measurements of the 7-ray emission will thus allow one 
to constrain the primary electron spectral index. 

4.3. VHE j-ray observations 

The Coma clust er has been observed i n the VHE 7-ray 
band with HES S (iDomainko et al.ll2007f) and with VER- 
ITAS (iPeridns] [20081). Both observation are consistent, 
and we focus on the HESS core observations, which pro- 
vide more stringent constraints on our model. The upper 
limit inferred from HESS observations on the flux above 
ITeV within 0.2° is F Coma < 8.3 • 10~ 13 phcm" 2 s" 1 . In 
fig. [7] we compare HESS's upper limit for the Coma clus- 
ter with the model predicted > 1 TeV 7-ray flux (within 
a disk of angular radius 9 as function of 9), assuming 
/inst^e = 0.01, /3 CO rc = 10~ 4 and B = 10fiG. Since the 
energy threshold for this measurement is very high, it 
provides robust constraints only on /3 C orc, through the 
predicted ir° decay luminosity. Only weak constraints 
on /inst??e may be obtained, due to the uncertainty in 
the maximal energy to which primary electrons may be 
accelerated, which is ~ few TeV. 

Using eqs. {27 ) . ([28" ]) and ([50 ]) . and 



0.2°/0.14° 



(1 



-3£c 



a + 1/2 



fdf ~ 0.43, 



we have 



oma, total ~6.8- 10" i5 /3 

core, —4 

+ 2.0- 10- 14 ((/ inst )o.2o?7 e )- 2 



ph 



(36) 



(37) 



(due to the strong dependence of the IC flux on the 
maximal energy of the primary electrons, we used the 
exact formulae for IC emission spectra, as described 
in § I2.5j) . The upper limit of HESS therefore implies 
/3 corc < 1.2 • 10 -2 (and (/ ms t)o.2°f7e < 0.42, with large 
uncertainty) . 



4.4. Radio observations 

In fig. [5] we compare the iThierbach et~al] £2003) com- 
pilation of flux density measurements of the Coma radio 
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Fig. 7. — The predicted > 1 TeV 7-ray flux within a disk of 
angular radius 6 as function of 8 for the Coma cluster, assuming 
/inst^e = 0.01, /3 C orc = 10 — 4 and B = lOttG. Line types are the 
same as i n fig. [T] The arrow re presents the upper limit reported 
by HESS (jDomainko et al.l l20071. which implies /3 COIC < 1.2 ■ 10 -2 . 
Only weak constraints on /i na t?7e may be obtained, due to the un- 
certainty in the maximal energy to which primary electrons may 
be accelerated, which is ~ few TeV. 
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Fig. 8. — A comparison of the [Tliicrbach etall (I2003T ) compi- 
lation of flux density measurements of the Coma radio halo with 
our model predictions, assuming B 3> BcmBi /Score = 2- 10 -4 , and 
that the radio halo flux is dominated by synchrotron emission from 
the secondaries (and taking into account the SZ effect). The spec- 
tral steepening observed above ~ 2 GHz is not robust (see text). 
Line types are the same as in fig. [T] 



halo with our model predictions, assuming B 3> -Bcmb, 
/3core = 2- 10 , and that the radio halo flux is dominated 
by synchrotron emission from the secondaries. We note 
that the spectral steepening observed above ~ 2 GHz is 
not robust, since it is the result of the substraction of 
two big numbers, the total flux and the flux of point 
sources, and the flux of point sources is not measured 
but rather extrapolated from lower frequencies assum- 
ing a constant spectral index. Since a steepening of the 
spectra is visible in other sources (for example, the two 
central galaxies of Coma: NGC4869 and NGC4874), by 
assuming a constant spectral index one overestimates the 
point source flux and u nderestimates the flux d ensity of 
the diffuse component ()Thierbach et al.l l2003f) . More- 
over, somewhat above ~ 3 GHz t he observed fl ux is sup- 
pressed also by the SZ effect (see lEnfilinll2002l for a de- 
tailed discussion). 
Using eq. (|34p and eq. (|2"9"|) we have 

f5„,coma — 5.7 • 10~ 1 /3 corc ,-4 Jy GHz. (38) 

Comparing with observations (see fig. [8]) this implies 
/3 coro ~ 2 • 10~ 4 . Note, that this constraint is consis- 
tent with the averag e value of /3 C orc — 10~ 4 derived by 
iKushnir et al.l (|2009T) for a complete sample of radio emit- 
ting galaxy clusters. 

5. ICM AND CR EVOLUTION 

In §|2] we derived the spectral and radial distribution of 
the nonthermal emission produced by ICM CRs for mas- 
sive clusters, M > 10 M©. We assumed that the frac- 
tion /3 CO rc of plasma energy carried by CR protons at the 
central regions of clusters, which dominate the emission 
from secondary particles, is nearly independent of clus- 
ter mass and that the scatter around its average value 
is small. We have argued, based on simple arguments 
and crude approximations, that /3 C oro — f/p/lOO. In what 
follows we present a simple model for the thermal his- 
tory of the ICM, including the effects of mergers, which 
allows us to obtain a more accurate estimate of the value 



of Pcr,p, its spatial dependence within clusters, and its 
scatter among different clusters. 

We first describe in § 15.11 our model for the accretion 
and merger history of clusters. In § 15.21 a simple model 
describing the effects of mergers is constructed. The re- 
sults of § 15.11 and § 15.21 are used in § 15.31 to construct a 
model for the evolution of CRs in the ICM. The effects 
of entropy changes that are not driven by gravity (i.e. 
not due to accretion and merger shocks) are discussed in 
§ 15.41 We assume throughout the analysis that the CR 
pressure is small compared to the therm al plasma pres- 
sure (as supported by observations, see IKushnir et al.l 
f2009h . 

5.1. Merger and accretion history 

We define cluster mass as the mass contained within a 
radius r2oo, within which the mean density is 200 times 
the critical density, M 2 oo,2 = M(r2oo) = (4/3)7rr2 00 x 
200p CI i t (z). We investigate the present, z = 0, non- 
thermal emission of M 200 ,o > 1O 14 5 M clusters. For 
each value of M2oo,0) chosen from a grid of masses 
within this rang e, we construct an ensem ble of " merger 
tr ees" using thelPress fc Schechterl (|1974h based scheme 
of lLacev k, Co le (1993). In order to incorporate the ac- 
cretion process into this scheme, we assume that merger 
events involving masses lower than a certain "resolution" , 
Mi, are accretion events. Since the mass accretion rate, 
M , calculated in t his m anner diverges within the scheme 
of lLacev fc Cold (|1993l ) as the time step tends to zero, 
we calculate M within the frame work of Press-Schechter 
theory (see § [B] for details) . Figure [9] shows M obtained 
in our scheme for different values of / acc = Mi /Afaoo,0) 
normalized to M2oo,z/tn.z, where tn. z is the Hubble time 
at redshift z. The accretion rate is averaged over 100 
realizations of a M 2 oo,o = 1O 15 M cluster. The figure 
illustrates that M depends weakly on the value of / acc , 
and is lower than M2oo. z /tH,z at low redshifts by a fac- 
tor of 2 — 5 for mass resolutions which are adequate for 
a description of a merger tree (i.e. that reproduce the 
general results of the Press-Schechter theory). The de- 
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Fig. 9. — Our estimate for the accretion rate, normalized to 
M200,z/tH,z, for different ratios / acc = M;/M 2 oo,0- The accre- 
tion rate is averaged over 100 realizations of a M2oo,0 = 1O 15 M0 
cluster. 



viation of M from M 2 qq,z /tH,z is expressed in terms of 
/i nst = MtH. z /M2oo,z- As discussed in § [3J 3D numeri- 
cal simulations indicate that the average value of /i ns t is 
« 0.5 at low z. We therefore choose M; = 2 • 10~ 3 M 2 oo,o, 
which yields /i nst ~ 0.5 for z < 1 and 0.5 < /i ns t < 1 
for 1 < z < 2.5. Our results depend only weakly on the 
value of /acc, and are little modified when / acc = 5 ■ 10~ 4 
is chosen instead of / acc = 2 • 10~ 3 . 

Next we describe how the gas profile is evolved with 
the addition of accreted mass. We assume that the halos 
are described, both b efore and after the a ccretion, by an 
NFW density profile (|Navarro et al.lll997l ). 



p(r) 



{r/r s )(l+r/r s y 



where p s — M s /(47rrf) and 



M, = 



M, 



J 20(1 



ln(l + r 2 Qo/r s ) - (r 2 oo/?'s)/(l + r 2 oa/r s ) 



(39) 



(40) 



r s is often expressed in terms of the "concentration pa- 
rameter", C200 = T2oo/r s . We adopt C200 = 4 for all of 
our systems, a value typical of galaxy clusters simulated 
in a ACDM concord ance cosmology. Numerous studies 
(e.g. lEke et ail |2001|) found that the mass dependence 
of the concentration parameter is weak, varying between 
different clusters by only a factor of ~ 2. 

We assume that the baryon density is proportional to 
the total (baryon and dark matter) density, p gas pa fop, 
where fo is the cosmic baryon fraction, in a greement with 
numerical simulations and observations (|Borgani et al.1 
l2004t iKravtsov et all [20051 : IVikhlinin et al.ll2006D . This 
assumption is not valid within cluster cores, where the 
gas density profile is flatter than that of the dark mat- 
ter. Furthermore, the NFW profile is different than the 
/3-model we used for the gas density in §[2J However, the 
choice of an NFW or a /3-model for the gas density is 
expected to have only a small effect on the inferred value 
of /3 CO re, due to the following reason. As we show below, 
the production of CR protons is dominated by accretion 



shocks and the evolution of their energy density is af- 
fected mainly by adiabatic compression. This, combined 
with the fact that both the NFW and the /3-model den- 
sity profiles imply a ratio of ~ 100 between the gas den- 
sity in the region that dominates the secondary emission 
and the gas density at ^200, imply that the modification 
of /3 CO rc due to the difference between the profiles of such 
models is small. 

Another inaccuracy introduced by using an NFW gas 
profile is an exaggerated cooling of CR protons in the 
very center of clusters due to the over estimate of gas 
density (see eq. ©). This, however, has only a minor 
effect on our results since the cooling affects only a small 
fraction of the gas mass (see § I5.3[) . We chose to work 
with an NFW density profile despite the above limita- 
tions, since it allows us to derive a very simple model for 
the description of cluster mergers (see § 15. 2p . 

In order to completely specify the properties of the gas, 
we must choose an entropy profile. Assuming hydrostatic 
equilibrium, 

the entropy is deduced from the equation of state, p = 

5 /3 

Kpgis ■ A boundary condition must be specified for 
eq. ([41 p . We choose to specify a value for the pressure of 
the gas at r 2 oo- Although there is some freedom in the 
choice of the pressure at r 200 , the entropy profile is rather 
insensiti ve to p(r = r 2 oo) (for p hysically reasonable val- 
ues, see iMcCarthv et aLll2007l) . Since at intermediate 
radii the NFW profile is well approximated by an isother- 
mal profile, i.e. p oc r~ 2 and d\nK/d\nM gas = 4/3, 
we choose a value of p(r 20 o) that establishes this near 
pure power-law e ntrop y distribution all the way out to 
r 200 . lVoit et all (|2002fl have demonstrated that groups 
and clusters are indeed expected to have near power-law 
entropy distributions out to large radii. 

Assuming that a hydrostatic NFW profile is preserved 
in the accretion process, we modify, following accretion, 
the density and the temperature of each mass element 
according to its Lagrangian location, in order to match 
the post-accretion hydrostatic NFW profile. We finally 
assume that the accreted mass was shocked at an infi- 
nite Mach number shock to the density and temperature 
at the outskirts of the post-accretion hydrostatic NFW 
profile. 

5.2. A simple model for mergers 

Let us next construct a simple model describing the 
evolution of ICM gas in merger events. We first note 
that the most frequent mergers occur between halos of 
similar mass. This is illustrated in fig. 1101 which shows 
the distribution of the ratio between the primary (heav- 
ier halo) mass, M p , and secondary (lighter halo) mass, 
M s , over 100 runs for M 2 oo,o = 1O 15 M . The distribu- 
tion is weighted by the secondary mass, since the ther- 
modynamic variables of the secondary halo are more af- 
fected in merger events than the primary ones. Accord- 
ing to fig. QUI high mass ratio mergers, M p /M s > 10, are 
rare. We therefore treat such mergers in a very approxi- 
mated manner, as described below. Note, that although 
the mass accumulated in high mass ratio mergers is a 
small fraction of the total cluster mass, it is not neces- 
sarily small compared to the cluster core mass, which 
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Fig. 10. — The distribution of primary to secondary halo mass ra- 
tio in mergers, obtained from 100 realizations of the merger history 
for M200.0 = 10 15 Mq. The distribution is weighted by the sec- 
ondary mass, since the thermodynamic variables of the secondary 
halo are more affected in merger events than the primary ones. 



dominates the secondary emission. Thus, high mass ra- 
tio mergers could have affected the nonthermal emission, 
had the mass of the secondary accreted halos been accu- 
mulated in the cluster core. Numerical simulations (e.g. 
iMcCarthv et al.l 120071 ) indicate that this is not the case: 
For M p /M s — 10 the secondary halo penetrates through 
the massive one and its mass is spread, on a Hubble time 
scale, over the entire primary halo. 

For the description of small mass ratio mergers 
we rely on detailed 3D numerical simulations (e.g. 
IMcCarthv et~aDl2007h . We assume that the halos are 
described, both before and after mergers, by an NFW 
density profile. We show that this assumption, combined 
with the inference from numerical simulations that dur- 
ing the merger the ICM is shocked twice (with the two 
shocks separated by adiabatic expansion) , completely de- 
termines the Mach numbers of the shocks and the magni- 
tude of the adiabatic expansion experienced by different 
fluid elements. 

We assume that the available energy to be thermal- 
ized in the merger is divided equally between the two 
shock episodes. This determines the Mach number of 
each mass element in each shock episode (for details, see 
IMcCarthv et al.ll2007l ). such that the only two free pa- 
rameters are the magnitudes of the adiabatic expansion 
experienced by each cluster. Assuming that the profile 
of the merged cluster is produced by sorting the gas el- 
ements of the initial clusters according to their density, 
the two adiabatic expansion factors are determined by 
minimizing the difference between the resulting profile 
and a hydrostatic NFW profile. In order to illustrate 
the procedure used, fig. [JTJ and fig. [12] show the resulting 
density and entropy profiles for a merger with a 1 : 1 
mass ratio. Merging the clusters, and sorting according 
to density leads to an exact NFW profile, with entropy 
which is lower than required by hydrostatic equilibrium. 
After the first shock, the entropy becomes higher but 
the resulting density exceeds the NFW density. Adia- 
batic expansion and a second shock then lead to density 
and entropy profiles which are very close to NFW at hy- 
drostatic equilibrium. The expansion factor required for 
obtaining the required final NFW density and entropy is 



Initial 

After 1st shock 

After adiabatic expansion 

After 2nd shock 



Q. 



0.2 



0.4 0.6 

M /M 

gas gas, tot 



0.8 



Fig. 11. — The gas density of the merged cluster as function of the 
accumulated gas mass for an equal mass merger, for a calibrated 
adiabatic expansion density decrease factor (see text) of ~ 0.17. 
The profile of the merged cluster is produced by sorting the gas 
elements of the initial clusters according to their density. 
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Fig. 12. — Same as fig. II II for the entropy of the merged cluster 



~ 6, corresponding to a density decrease factor of ~ 0.17 
It is remarkable that this simple model reproduces a 
post-merger hydrostatic NFW halo. We note that assum- 
ing a different density profile for the halo (for example, a 
(3- model, see § d|) does not allow one to preserve the den- 
sity profile in the post merger cluster using such a simple 
model. The density decrease factors determined by this 
method for mergers of different mass ratios are shown in 
fig. rj5] Above a mass ratio of 10 : 1 we simply sort the 
gas elements according to their density, and change the 
density and the temperature of each mass element ac- 
cording to its Lagrangian location, in order to reproduce 
the post-merger hydrostatic NFW profile. The details of 
large mass ratio mergers have only a small effect on our 
results. 

One of the results of our model is the merger shocks' 
Mach number distribution. We show this distribution, 
weighted by the mass of the shocked gas, in fig. [21 
based on 100 merger history realizations for a cluster 
with Af2oo, o = 1O 15 M (the first bin is not shown, since 
almost all of the shocks are included in it). As can be seen 
from the figure, the maximal Mach number achieved in 
merger shocks is smaller t han 2.5, in agre e ment with the 
numerical simulations of Skill man et al.l (|2008l) . Since 
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Fig. 13. — Calibrated adiabatic expansion factors as function of 
merging clusters mass ratio. 
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Fig. 14. — Merger shocks' Mach number distribution, weighted 
by the mass of the shocked gas, averaged over 100 merger history 
realizations, for a cluster with M2oo,0 = 10 15 Mq. The first bin is 
not shown, since almost all of the shocks are included in it. The 
maximal Mach number achieved in mergers is smaller than 2.5. 



such weak shocks are assumed to produce very steep CR 
spectra (see § Q}, with normalization that depends on 
poorly known conditions at the injection energy of the 
CRs, we ignore acceleration in merger shocks hereafter. 

5.3. CR evolution 

After determining, in § 15.11 and § 15.21 the thermody- 
namic history of the ICM plasma, we turn now to the 
evolution of the CR population. We follow here only the 
evolution of CR protons, since we are only interested in 
the value of (3 COTC . The processes that we include in the 
calculations are: injection of CRs in accretion shocks (we 
neglect injection in weak merger shocks, see § 15. 2p . proton 
energy losses due to Coulomb and i nelastic nuclear colli - 
sions, using the parametrization of iKamae et al.l (2006) 
for the p-p interaction, modification of the CR energy 
density due to adiabatic expansion/compression (assum- 
ing the CRs to behave as a relativistic gas with an adi- 
abatic index of 4/3). As explained in the introduction, 
we assume that CR diffusion is not important, and that 
the CRs are coupled to the thermal plasma by magnetic 
fields. We also assume that the CR pressure is small 
compared to the the rmal plasma pressur e (as supported 
by observations, see lKushnir et all 12009). We calculate 
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Fig. 15.— Pcnhv at 100 GeV (blue) and at 10 GeV (black), as 
function of the accumulated mass for a typical Af2oo,o = 1O 15 M0 

cluster. The red line shows a p gas scaling of /^cr/^pi normalized 
at the accretion shock. 



numerically the CR content of each mass element in the 
cluster, following the thermodynamic evolution of this 
mass element and describing the CR energy distribution 
using a discrete (Lagrangian, logarithmically spaced) set 
of energy bins. 

Figure [15] presents typical distributions of /3cr,p(s = 
100GeV)/?7p and /?cr,p(£ = 10GeV)/r? p , obtained for 
one realization of the evolution of a 10 15 Af© cluster. 
/?cr,p depends weakly on energy, as demonstrated in 
the figure, from ~ 10 GeV up to the maximal energy 
to which the CR protons are accelerated. The spatial 

— 1/3 

distribution of /?cr,p approximately follows /9 gas , with 
normalization fixed at the accretion shock. /?cr,p falls 
somewhat below this scaling at the cluster center mainly 
due to the effect of weak shocks (which lead to compres- 
sion without significant production of high energy CRs) 
and due to proton energy loss at the high density cluster 
core. The latter effect is exaggerated in our calculations 
due to our assumption that the the gas density follows 
the dark matter density. However, cooling affects only a 
small fraction (1%) of the gas mass. The fact that Pcr, p 

— 1/3 

follows approximately p ga J implies that its evolution is 
determined mainly by the adiabatic compression of the 
ICM plasma, and that mergers do not significantly mod- 
ify /?CR,p- 

We define (3 COIC as the mass average of /3ck at 100 GeV 
within the inner 10% of the cluster mass, which cor- 
responds to the core region. Figure [16] shows the av- 
erage value and the scatter of f3 COIC /r]p as function of 
the cluster mass at z = 0. As can be seen from the 
figure, Pcore/Vp — 1/200 is a good approximation for 
-^200,0 > 1O 14 5 M . (3 COIC depends only weakly on clus- 
ter mass, changing by a factor of less than 2 over one 
decade of cluster mass, with a factor ~ 2 scatter between 
different clusters of given mass. 

5.4. The effects of "non- gravitational" entropy changes 

Our model for the ICM evolution takes into account 
entropy changes of the ICM driven by gravity only (ac- 
cretion and merger shocks), and assumes that the galaxy 
clusters' density and entropy profiles are self-similar. It is 
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Fig. 16. — /3core (the mass average of /3cr, p at 100 GeV within the 
inner 10% of the cluster mass) as function of cluster mass at z = 0. 
The average value (for a given cluster mass) was obtained from 100 
realizations of merger histories, and the error bar represents the la- 
scatter of the Pcn,p distribution. 



well known (see e.g. lArnaud &: Evrard1l999( ) that models 
taking into account only gravitational effects and assum- 
ing self-similarity lead to a correlation Lx oc T 2 between 
the bolometric X-ray luminosity, Lx, and the tempera- 
ture of the clusters, T, which is fla tter than the observed 
corre lation, Lx oc T 2 5 (see, e.g. iReiprich fc Bohringerl 
2002). This deviation from self-similarity is commonly 
explained as due to entropy increase of the gas at 
some high redshift by non-gravitational processes such 
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1999). Although such entropy increase does not modify 



significantly the entropy profile of the massive clusters in 
which we are interested, one may worry about its pos- 
sible effect on the CR population, since part of the gas 
residing in high mass clusters has, in the past, been as- 
sociated with low mass clusters. The CR energy density 
in low mass clusters may be reduced by an early entropy 
increase, which reduces the Mach numbers of accretion 
shocks. Repeating our calculations including an early en- 
tropy increase we find, however, that this effect is rather 
small. Early (non-gravitational) entropy increase may 
reduce the Mach numbers of accretion shocks to values 
which significantly affect CR production only for clus- 
ters of mass < 10 13 M Q , and the suppression of CR pro- 
duction in such shocks reduces the the value of /3 CO ro in 
massive clusters by no more than 20%. 

Another non-gravitational effect that may modify the 
ICM entropy is radiative cooling at the cores of massive 
clusters. This effect may lead to an increase of the ratio 
of CR to thermal ICM energy density and may therefore 
affect the predicted secondary emission. However, nu- 
merical simulations (|Pfrommer et al.l 120071 ) predict that 
the effect of cooling on the nonthermal secondary emis- 
sion is not very large, as can be estimated by examining 
the simulated 7-ray emission from massive clusters, for 
radiative and non-radiative simulations. This is proba- 
bly due to the fact that cooling affects mainly the inner 
few tens of kpc, while the secondary emission is dom- 
inated by a larger region, ~ 300 kpc. We note that it 
is well-known that the effects of cooling flows are exag- 



gerated in numerical simu lations fe.g. lBalogh et alj|2001t 
iBorgani et al.i r2004. 200(| and references therein). 

6. DISCUSSION 

We have derived analytic expressions, eqs. ([9jl-([26]). 
which approximately describe the spectral and radial dis- 
tribution of the nonthermal emission produced by (pri- 
mary and secondary) CRs in massive, M > 10 14 ' 5 M Q , 
galaxy clusters. These expressions depend on the (ob- 
served) cluster thermal X-ray properties and on two 
model parameter, f3 cmc and r\ e . We have shown (see 
§ [5]) that CR production is dominated by strong accre- 
tion shocks, and that the energy density of CR pro- 
tons is mainly affected (after production) by the adi- 
abatic compression of the ICM plasma. Simple argu- 
ments based on crude approximations imply that adia- 
batic c ompression shou l d lead to /3 CO re ~ r^/lOO (jj [2j 
see als dPfrommer et all 120071: iJubelgas et al.ll2008f ). A 
more detailed analysis of the thermodynamic history of 
the ICM (§ [5]) yields /3 CO re — %/200, nearly indepen- 
dent of cluster mass and with a scatter Aln/3 corc ~ 1 
between clusters of given mass (see fig. IT6]) . The analy- 
sis presented in § [3] includes a simple description of the 
thermodynamic history of the ICM, as determined by 
merger and accretion events, and a calculation of the en- 
ergy and spatial distribution of CR protons. We have 
shown in § [4] that our model's results agree with those 
of detailed numerical calculations, and that discrepan- 
cies between the results of various numerical simulations 
(and between such results and our model) are due to 
inaccuracies in the numerical calculations. 

The charged secondary emission depends on the the 
strength of the magnetic field in the cluster core. It 
is given in § 12.31 in terms of B/Bqmb, the ratio of the 
magnetic field to Bqmb, defined as the magnetic field 
for which the magnetic energy density equals the CMB 
energy density, £?cmb — 3/iG. Since the secondary 
emission is dominated at high energy, > 1 eV, by sec- 
ondary 7T° decay and by IC emission of primary elec- 
trons, the uncertainty in the value of B affects only the 
predicted radio emission, which is not the main focus of 
the current paper. The d etailed discussion of r adio emis- 
sion from clusters given in lKushnir et al.l (j2009T ) indicates 
that B is within the range of ~ 1/iG to ~ 10/iG. 

Our model predicts that the HXR and 7-ray lu- 
minosities produced by IC scattering of CMB pho- 
tons by electrons accelerated in accretion shocks ex- 
ceed the luminosities produced by secondary parti- 
cles by factors ~ 500(?7 e /77p)(r/10keV)~ 1 / 2 and ~ 
150(ry e /r? p )(T/10keV) _1 / 2 respectively, where T is the 
cluster temperature. Secondary particle emission may 
dominate at the radio and VHE (> 1 TeV) 7-ray bands. 
Our model predicts, in contrast with some earlier work 
that neglected the primary IC emission, that the HXR 
and 7-ray emission from clusters of galaxies are extended, 
since the emission is dominated at these energies by pri- 
mary (rather than by secondary) electrons. Our model is 
supported by observations of the Coma cluster, where the 
HXR image does not show any significant substructure, 
the measured HXR flux corresponds to a reasonable ef- 
ficiency of electron acceleration, rj e ~ a few percent, and 
the radio emission corresponds to a reasonable efficiency 
of proton acceleration, f3 COIC ~ 10~ 4 implying r\ v ~ a 
few percent (see § |4j) . As explained in § 14.11 (see also 
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iKushnir et al.ll2009l) . in order for secondary emission to 
reproduce the HXR emission of Coma, an unreasonably 
high CR energy density is required. 

We have assumed in our analysis that the ICM plasma 
is isothermal and in hydrostatic equilibrium. Devia- 
tions from this simple model near the virial radius may 
change our estimates for r\ e . p . Although such deviations 
are only wea kly constrained by obs ervations, both obser- 
vational (e.g. iVikhlinin et al.ll2005| ) and theoretical (e.g. 
iRoncarelli et al.ll2006f) analyses indicated that they are 
not large (for example, the accretion shock temperature 
is lower than the virial temperature by no more than a 
factor ~ 2). Modifications of the ICM properties near the 
virial radius may be easily incorporated into our model. 
Improved (observational) determination of the ICM pro- 
file near the virial radius will therefore allow one to im- 
prove the accuracy of the determination of We,v 

Our model predicts, that future HXR observations 
(e.g. NuStar, Simbol-X) and space-based 7-ray observa- 
tions (Fermi) will lead to detection of clusters of galax- 
ies as extended sources. Since EGRET's sensitivity is 
marginally too low for such detection (see §E]), we pre- 
dict Fermi's ~ 50 times better sensitivity would suffice. 
Single sources could be detected and resolved in both en- 
ergy bands (the angular resolution of Fermi can be as low 
as 0.1° and that of future HXR observation may reach 
tens of arcsec) . We give in table Q] the predicted flux 
above 50 GeV within 0.1° for the 10 most luminous clus- 
ters fr om the extended sample of iReiprich k, Bohringerl 
(2002), assuming r] e fi nst = 0.01. We also give the cluster 
angular radius. For cluster diameters larger than 0.2° it 
would be possible to distinguish between extended and 
core emission with Fermi, thus discriminating between 
secondary emission from the core and primary emission 
from the accretion shock. 

Figure [7] shows that the current upper limit on the 
VHE 7-ray flux from the core of the Coma cluster is 
higher than the predicted flux by a factor of ~ 100. This 
implies that a detection of the predicted Coma core flux 
is not possible with current generation imaging Cerenkov 
telescopes. Imaging telescopes may, however, allow us 
to detect the predicted high energy emission produced 
by the accretion shock. It should be pointed out that 
lowering the threshold energy of the imaging Cerenkov 
telescopes to ~ 0.1 TeV is important for this type of 
observation, since there is a large uncertainty in the pre- 
dicted flux at energies > TeV, due to the uncertainty in 
the exact value of the maximal energy to which electrons 



can be accelerated, which is ~ few TeV. The sources in 
table [T] are also optimal for imaging telescopes. 

In order to detect pion decay emission from cluster 
cores, two approaches may be adopted. One may con- 
struct high resolution (~ 0.1°) HXR maps and compare 
them to the 7-ray maps. Since the correlation between 
the HXR and the 7-ray maps should be due to the pri- 
mary IC emission from the accretion shock, their differ- 
ence should trace the pion decay contribution. Alterna- 
tively, one may identify clusters which may be resolved 
and for which the surface brightness of the core is dom- 
inated by pion decays, so that a sharp decrease in the 
surface brightness with radius may be observed. Tabled] 
gives a list of the best candidates for such detection: The 
predicted flux above 50 GeV within 0.1° is given for the 
10 clusters with t he highest pion to IC ratio from the 
extended sample of lReiprich fc Bohringerl ()2002f ) , assum- 
ing T] e finst = 0.01 and /3 C orc = 10~ 4 . We considered only 
clusters with angular radius exceeding 0.2°, which may 
be resolved. 

As illustrated by the analysis of the observations of 
the Coma cluster, high energy and radio observations 
of a controlled sample of clusters will allow one to cali- 
brate the model parameters ry e /i nst and /3 C oro- Determi- 
nation of /i ns t and Pcarc/Vp from numerical simulations 
would then allow one to determine r\ v and rj e . Analy- 
sis of cluster radio emission, leading to a determination 
of /3 CO ro and of the stre ngth of magne t ic fiel ds in clus- 
ter cores, is reported bv IKushnir et al.1 (|200 9). Analysis 
of all clusters observ ed in HXRs, aimed a t deter mining 
We f inst, is reported bv IKushnir fc Waxmanl (|2010f ). Since 
the group of clusters observed in HXRs is not a complete 
controlled sample, the determination of Tfe/inst based on 
current observations may be biased. A controlled sam- 
ple may be produced by future HXR observations (e.g. 
NuStar, Simbol-X). Moreover, such missions may pro- 
duce high resolution HXR maps of clusters, testing the 
prediction that the emission is extended, i.e. originating 
from the accretion shock. Finally, we note that measure- 
ments of the nonthermal emission in different bands (e.g. 
HXR and 7-rays for electrons, radio and VHE 7-rays for 
protons) would allow one to constrain the energy distri- 
bution of the accelerated particles. 
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APPENDIX 

A. EMISSION MECHANISMS 
Thermal emission 
The bolometric brcmsstrahlung emissivity at some radius is given by 

e x (r) * ^o-TaecV^T^m-^e 1 + ^) f^(r), (Al) 

(assuming the thermal Gaunt factor to be 1 and neglecting elements heavier than Helium, which may increase the 
bremsstrahlung emissivity by a few tens of percent). Here a e is the fine structure constant and we use the definitions 
/(M 

/(fiHTUp) and n Hc = P g a S /(m m p ) for the electron, Hydrogen and Helium number 
densities, respectively (which yield fi e ~ 1.14, hh — 1-33 and ^tHe — 16 for fully ionized gas with hydrogen mass 
fraction x == 0.75). 
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Table 1. The predicted IC flux above 50GeV within 0.1° for. 

THE 10 MOST LUMINOUS CL USTERS FROM THE EXTENDED SAMPLE OF 
R.EIPRICH <fc BOIIRINGER (2002), ASSUMING rj e f inst = 0.01. 



Cluster name Flux [10 12 phcm 2 s 1 ] Angular radius [deg] 



A2163 


0.71 


0.21 


A1914 


0.45 


0.22 


Ophiuchus 


0.40 


1.3 


A3888 


0.38 


0.25 


A0754 


0.34 


0.65 


A1689 


0.33 


0.18 


Triangulum 


0.31 


0.63 


A2142 


0.31 


0.35 


A3266 


0.27 


0.56 


A2029 


0.27 


0.40 



Table 2. The predicted flux above 50GeV within 0.1° for the 10 clus- 
ters with the highest pion t o ic ratio from the extended sample of 

ReIPRICH & BOHRINGER (2002), ASSUMING rj e f inst = 0.01 AND /3 C ore = 10~ 4 . 



Cluster name Pion decays flux IC flux ratio Angular radius 





[10- 13 phcm- 2 s- 1 ] 


[10- 12 phcm- 2 s- 1 ] 




[deg] 


Perseus 


3.1 


0.15 


2.0 


1.4 


A3526 


0.49 


0.05 


1.0 


1.7 


NGC 4636 


0.026 


0.003 


0.80 


2.2 


2A 0335 


0.30 


0.038 


0.79 


0.50 


NGC 5813 


0.017 


0.0024 


0.73 


1.3 


Ophiuchus 


2.7 


0.4 


0.67 


1.3 


A0262 


0.12 


0.018 


0.65 


0.80 


A2199 


0.41 


0.072 


0.56 


0.72 


A0496 


0.33 


0.059 


0.56 


0.57 


NGC 5044 


0.034 


0.0061 


0.56 


1.1 



Neutral secondaries emission 

The p-p 7-ray emissivity per logarithmic photon energy bin, ve pp , due to the decay of ne utral pions produced in 
inelastic nuclear collisions, is given to a very good approximation by ()Katz fc Wa xman 2008) 

1 



:2f m> o.i^ ( — + — ) 

\cfey m p \fi H MHc/ 

cpcnm- 2 /,- 1 (— + —) Tp 2 glis (r)a™\ (A2) 



where cr p n p el ~ 40 mb and / pp ~ 0.75. Since the CR protons have flat spectra, ve pp is not energy dependent at the 
relevant photon energies, from ~ 0.1 GeV (eth — 1-22 GeV is the threshold energy for pion production and the photon 
energy is ~ 0.1 of the proton energy) up to the cutoff energy due to pair production in in teraction with the IR 
background, which should be around ~ lOTeV for nearby clusters (|Franceschini et aL1l2008l) . Using eq. (|A1[) and 
eq. (|A2[) . the total p-p 7-ray luminosity, vLf, p , is related to the bolometric bremsstrahlung luminosity, Lx, by 

1 _l_ 1 3 n 9 f /r^ nc ^ R 

VK 1 4 ~7= 1 L X 

" Mir A»Ho \ -rzCTTOt-e^/ m e C 2 



37T k 

,11 , rj.1/2 ( Lx \ I 



Using the correlations eq. (JSJ), we have 

vLl? -1.2- 10 41 /3 core ,- 4 Ii 3 06 ergs- 1 . (A4) 
The p-p 7-ray surface brightness above some energy £„ !Bl i n a t some distance f = r/r c from the cluster center is given 
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by (for > 0.5, see lSarazin fc BahcaUlfl97l 

^P, min (r)^^0.886^LPPr c 7 2 (3/? - (l + f 2 ) ^ +1/2 (max( e „, min , OA^))' 1 



max(e l/jmin ,0.1eth) V d , s 2\ -3/3+1/2 _ 2 -l -l / » r\ 

' (1+r ) phcm s sr (A5) 



lOGeV 

(neglecting the 77 -+ e + e~ cutoff). Using the correlations eq. (JB]) we have 

^ min (f) ~ 4.1 • 10- 7 (3/3 - |) /W^ ' 42 ( ^^GeV^O ' (! + r 2 )~ 30+1/2 phcm- 2 s- 1 sr- 1 . (A6) 

/C emission from charged secondaries 

p-p collisions also produce secondary electrons and positrons, which cool by emitting synchrotron radiation and by 
IC scattering of CMB photons. We assume that the distribution of the secondaries is in a steady state, where in the 
relevant energy bands the secondaries that are generated lose all their energy to radiation. Since the CR protons have 

flat spectra, the secondary production (eL\ ) in the range O.leth < £ e ± < 0.1e m ax is not energy dependent. To a very 
good approximation, the production in this energy range satisfies 

sLf ~ /e ~ ± /e+ vL*?, (A7) 

where f e - ~ 0.8 and / e + ~ 1.2. As explained in § we consider only secondaries with cooling time shorter than 
idyn- For values 1 — 10 /iG of the magnetic field, the Lorentz factor of secondaries with cooling time which equals the 
dynamical time is 200 < 7 CO oi S 2000. 

To a very good approximation, in the range of photon energies 7 2 oo1 3Tcmb(1 + -z) 4 < e p h < 7^ ax 3TcMB(l + z) 4 , 
where 7 max — 0.1e max /m e c 2 (e ma x is the maximal energy of the CR protons and the secondaries energy is ~ 0.1 of the 
proton energy), the luminosity due to IC scattering from the secondaries is given by 

j^IC.e* ^ 1 c ^e ± ^CMB ^ fe+ + fe~ v j^VP ^CMB 



and the surface brightness is given by 



2tt 5 /2 • " c V 2 



^ ^CMB a j.1/2 / £ 

b cmb + B 2Pcmc '- 4 1 Uf 2 3 • lO^ergs-V Ura 200 k P c 

/ £V,min \ 4 /-. _2\— 3/9 + 1/2 _ 2 —1 —1 

x ( 10GcV J ( 1 + r ) P hcm s sr ' ( A9 ) 

Using the correlations eq. (|6]), we have 

vL^ ~ 3.0 • 10 40 „^mb /3corc.-4r i 3 - 06 erg s' 1 , (A10) 

CMB ' 

and 

Sl C ^tjf) ~ 1.0 • 10- 7 (3/3 - |) ^l^&^Tr 2 (^) _1 (1 + f 2 )^ 3 ^ +1/2 phcm- 2 s^ 1 sr" 1 . (All) 



Synchrotron emission from charged secondaries 

To a very good approximation, in the range of photon energies 7 2 00 i£o < £ph < 7max e 0i where £o = 3heB / '47rm e c, 
the synchrotron luminosity from the secondaries is given by 



R 2 

CMB 



3-3 -10 40 a B \ R cmc ^Tl' 2 ( L f r ) (1 + z)- 4 ergs- 1 . (A12) 

b cmb+ b \h 70 3-10 45 ergs l J 
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and the surface brightness is given by: 

sr c < e >)-^o.886,Lr c < e V (3/3- 1) (l+fT^V 1 

3\_B^ a / Lx 



2/ B% MB +B 2 ^ ie ' 1 Wo 2 3- 10 45 ergs-V V^q 1 200kpc, 

x f ^— Y 1 (l + f 2 ) _3,3+1/2 (l + z)- 4 mJyarcmin- 2 . (A13) 

V 1.4 GHz/ 

Using the correlations eq. (|6]), we have 

^sync.ei _ 3 q . 1Q 40 bPq^ ^^{l + z)' A ergs" 1 (A14) 



J CMB 

and 



r c,e± (r) ^ 9.8 ( 3/3 - J ) R2 " , R2 Pco ie ,-AT?- i2 

\ 1 J CMB ' 

x (_^_) _1 (l + f 2 )" 3 ' 3+1/2 (l + z)- 4 mJyarcmin- 2 . (A15) 

Primary electrons IC emission 

As in the case of the charged secondaries, we assume that the distribution of the primaries is in a steady state, since 
electrons (at the relevant energies) lose all their energy to radiation on a time scale short compared to the cluster 
dynamical time, ta yn ~ 1 Gyr. Unlike the secondary electrons and positrons, which lose energy through both IC and 
synchrotron emission, primary electrons lose their energy mainly by IC scattering of CMB pho tons, since the magnetic 
field at the accretion shock is expected to be weak, ~ 0.1 fiG <C -Bcmb ()Waxman fc Lo cb 2000). Thus the steady state 
assumption holds for primaries with Lorenz factors 7 > 7 CO oi ~ 2000. 

To a very good approximation, in the range of photon energies 7 2 oo1 3Xcmb(1 + z) 4 < e p h < 7^ ax 3TcMB(l + z) A , the 
luminosity of IC emission by primary, shock accelerated, electrons is given by 

^ic.shock _ 1 3 Vehiump)-^ ^ Mgoo 
2 2 ln(p max ) - ln(m e c) m " t H 

~1.8-10 43 (/ inst7 7e)_ 2 ^ 3/2 (J^\T^ /2 Z(z)evg S -\ (A16) 

where (/inst'?e)_2 — /inst?7e/10 _2 , fb = ^b/^m and Z[z) = (tffH(z)) . Assuming that the primary electron emission 
originates from a thin layer of thickness w behind the accretion shock, the IC surface brightness is given by 



3\ B 



S%£?(r) - 3~ 2~ 2 — £ ir/r 200 ,w/r 20 o) ^ c > sh ° ck (e^y 1 



1.5-10- 7 (f inst rj e )_ 2 ^ 2 (^ 17 
) £,{ r / r 200,w/r 200 ) 



b_ \ j,3/2 



.lOGeV/ 

x Z(z)^ (z)phcm" 2 s~ 1 sr- 1 , (A17) 



where 

3(vT^2-- v /(l-y) 2 -ic 2 ) 



a*,y) = < , i-u-») a ~' x ~ 1 y ■ (Ai8) 

l-y<x<l 

For w <C ^200; y) can be approximated in the regime x < 1 — y as 

ax,y)c±^=L=. (A19) 
Vl — x z 

The thickness of the emitting region is approximately given by the product of the cooling time of the emitting electrons 
and the velocity of the downstream fluid relative to the shock velocity. As explained in §[21 the approximation w -C r 2 oa 
holds for 7 > 7 cool — 2000. 
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Primary electrons synchrotron emission 
The synchrotron luminosity from these electrons in the range of photon energies 7 2 ool £o < £ph < 7max £ 0j 1S given by 

B 2 



vL s J 



nc. shock , 



vL 



I C, shock 



f>2 

CMB 



~1.7-10 40 (/inst?7e)_ 2 /? 3/2 



fb 



0.17 



T 1 5/2 B2 7 Z(z)(l + z)- 4 ergs- 



(A20) 



where £>_7 = B/0.1 fiG. The synchrotron surface brightness from the shock is 



gsync 



shock 



(r): 



8irh 



< (r/r 200 , w/r 200 ) vL s J™' shock (z/)" 1 



200 



~1.4- 10- 2 (/instr?e)_ 2 ^ 1/2 



A 



0.17 



Tl ,2 Bi 7 



1.4 GHz 



x £ ( r / r 200, w / r 20o) z (z)h 2 (z)(l + z) 4 mJyarcmhT 
B. PRESS-SCHECHTER ACCRETION RATE 



(A21) 



In the Press-Schechter formalism, it is convenient to use u> = S c (z) = 5 c q/D(z) as the time variable (where <5 Cj o is 
the density contrast estimated by linear theory for the virialization of a spherical halo, and D(z) is the linear growth 
factor), and S(M) = a 2 (M), the mass variance, as the mass variable. The probability for a step AS in a time-step 
Aw is 



P(AS, Au)dAS 



1 



Aw 



2^ (AS) 3 /- 



■ exp 



(Au^y 

2AS 



dAS. 



(Bl) 



Changing variables, x = Aw/(2v / A5), this becomes a Gaussian distribution in x with zero mean and unit variance. 
Since AS oc Aw 2 , the accretion rate, AS 1 / Aw, is linear in Aw and thus tends to zero as Aw — > 0. 

In order to obtain a finite mass accretion rate in our scheme, we introduce a slight correction to the: treatment of 
masses below the mass resolution, Mf. Each time a mass step AM is drawn, which is below the mass resolution Mi, we 
assume that mass was accreted during Aw at a constant rate, (dM / doj) Mi acc , which is the average accretion rate below 



Mi. Let us next calculate (dM/du 



'M;,acc- 



The average mass of the progenitors of some mass Mq that contributed to 



M through accretion over a time step Aw is 

(■M 



P(M\M , Au)MdM = M 



Since 



we have 



M a -Mi 



1 - erf I 



Aw 



^2S{M l ) - 25(M ) 



Mo 



Aw a^o V 7T y/S{M t ) - S(M ) 



fdM_\ 

V duj / M( ,acc 



7T s jS(Mi) - S(M Q ) ' 



(B2) 



(B3) 



(B4) 
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